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Abstract 


Spurred by rich, multi-wavelength observations and enabled by new 
simulations, ranging from cosmological to sub-pc scales, the last decade 
has seen major theoretical progress in our understanding of the circum- 
galactic medium. We review key physical processes in the CGM. Our 
conclusions include: 

W The properties of the CGM depend on a competition between 
gravity-driven infall and gas cooling. When cooling is slow 
relative to free fall, the gas is hot (roughly virial temperature) 
whereas the gas is cold (T ~ 10* K) when cooling is rapid. 

E Gas inflows and outflows play crucial roles, as does the cos- 
mological environment. Large-scale structure collimates cold 
streams and provides angular momentum. Satellite galaxies 
contribute to the CGM through winds and gas stripping. 

W [n multiphase gas, the hot and cold phases continuously ex- 
change mass, energy and momentum. The interaction between 
turbulent mixing and radiative cooling is critical. A broad 
spectrum of cold gas structures, going down to sub-pc scales, 
arises from fragmentation, coagulation, and condensation onto 
gas clouds. 

W Magnetic fields, thermal conduction and cosmic rays can sub- 
stantially modify how the cold and hot phases interact, al- 
though microphysical uncertainties are presently large. 


Key open questions for future work include the mutual interplay be- 
tween small-scale structure and large-scale dynamics, and how the 
CGM affects the evolution of galaxies. 
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1. INTRODUCTION 


Observations indicate that galaxies from dwarfs to massive ellipticals are enclosed by mas- 
sive gaseous atmospheres, known as the circumgalactic medium (CGM). The total gas mass 
and the metal mass in the CGM can exceed the corresponding masses in galaxies (e.g. [Tum]] 
[linson et al.2017). Moreover, it has become clear in recent years that the CGM crucially 
affects the evolution of galaxies by mediating interactions between galaxies and the larger- 
scale intergalactic medium (IGM). Gas inflows from the cosmic web are necessary to sustain 
star formation in galaxies over cosmological timescales, while galactic winds play a critical 
role in regulating star formation rates (SFRs). Studies of the CGM therefore constrain, 
or make predictions for, the mass distribution, kinematics, thermodynamics, and chemical 
abundances of the gas flows that regulate galaxy formation. 

An earlier Annual Reviews article by [Tumlinson et al.| (2017) concluded that the CGM 
contents are now reasonably well characterized observationally and that key questions going 
forward include the physics that govern the CGM and how it interacts with galaxies. In 
the last decade, there have been a number of theoretical developments directly relevant to 
answering these physics questions. On large scales, the major advances include cosmological 
hydrodynamic simulations that now produce broadly realistic galaxy populations and which 
have been used to analyze CGM gas flows on large scales (for a recent review of cosmological 
simulations of galaxy formation, see E 2020). On small scales, there 
has been similarly important progress studying processes that are not well resolved in 
cosmological models, including the microphysics of how cold and hot gas phases exchange 
mass, momentum, and energy, and the inclusion of physics beyond ideal hydrodynamics, 
such as magnetic fields, thermal conduction, and/or cosmic rays. 

Our goal in this article is to draw on these recent advances and summarize our current 
understanding of the key physical processes that operate in the CGM. Our point of view is 
primarily theoretical, but our choices of topics are in many instances motivated by obser- 
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Figure 1: Schematic illustration of important properties of the CGM. The CGM is a complex 
system involving physics on an enormous dynamic range, from the dark matter halo on scales 
2, 100 kpc to sub-parsec structure in cold gas. Key concepts include the interplay between 
gas cooling, heating, and gravity in the halo and how the hot and cold phases exchange mass, 
energy, and momentum. The bottom panel shows a central galaxy whose star formation 
is fueled by a mixture of cold (blue) and hot (yellow) accreting gas and which powers a 
multiphase galactic wind, while the top two panels zoom onto a highly structured cold cloud 
complex (left) and a turbulent mixing layer (right). A complete understanding of the CGM 
and accurate observational predictions require consideration of various additional processes, 
including the effects of satellite galaxies (as in the bottom right) as well as magnetic fields, 
thermal conduction, cosmic rays, and feedback from accreting black holes (not shown here). 
Figure credit: Aaron M. Geller (Northwestern University/CIERA/IT Research Computing 
Services). ©2022 Aaron M. Geller. 
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vations. An understanding of CGM physical processes is relevant to several key questions, 
including: How does gas flow into CGM and accrete onto galaxies? How do galactic winds 
affect the CGM? How does the CGM affect the formation and evolution of galaxies? How 
is the multiphase structure of the CGM, which is critical to the interpretation of many 
observations, produced? What are the important physical scales in the CGM and what 
are the requirements to produce realistic simulations? We envision that our audience could 
range from new graduate students entering the field to more experienced CGM researchers 
interested in a summary of recent theoretical developments. 

By CGM, we typically refer to the gas within one virial radius Rvir of dark matter halos, 
but outside galaxies. We stress, however, that CGM processes such as galactic outflows can 
reach larger radii and we do not exclude such gas only because it has crossed the somewhat 
arbitrary virial-radius boundary. We focus mainly on the CGM around isolated galaxies, 
in halos of total mass up to a few times 10'? Mo, corresponding to central galaxies of 
order the mass of the Milky Way (the L, mass scale), or a few times this value. The most 
massive dark matter halos (My, > 10? Mo) correspond to galaxy groups and clusters of 
galaxies and they host the intra-group medium (IGrM) and the intracluster medium (ICM), 
respectively. While some physical processes are common to the CGM and the IGrM/ICM, 
we generally avoid discussing processes that are specific to group or cluster environments. 


There are other review articles which focus on groups and clusters (e.g.|Kravtsov & Borgani 
2012| |Donahue & Voit|2022). 


We organize our review into two main parts: one on cosmological processes (9) and one 
on small-scale processes (§3). The main theme of the section on cosmological processes are 
the properties and physics of gas flows on the scale of the halo, which connect the IGM to 
galaxies. The main theme of the section on small-scale processes is physical processes that 
arise in a multiphase medium, particularly how gas flows between different phases. A key 
emphasis is on processes that drive the formation, destruction, and structure of cold gas. 
The properties of this cold gas are important because common observational techniques, 
including as rest-UV absorption and emission, are sensitive to the cold gas phase. Section [4] 
combines cosmological and small scales in a discussion of the requirements for resolving cold 
gas and of the prospects for modeling cold gas in cosmological simulations. We summarize 
our outlook and outline key areas for future research in 

The artist’s conception of the CGM in F igure [I] previews some key themes. The figure 
illustrates a huge dynamic range and the importance of hot and cold phases. This illustra- 
tion contains several of the same concepts as [Tumlinson et al) (2017)’s Figure 1, which has 
been widely used to summarize key CGM processes, including filamentary accretion and 
outflows. We have produced a new cartoon picture to emphasize some of the complexity 
expected in the CGM, including messy structure on scales ranging from the halo (= 100 
kpc) to turbulent mixing layers (< 1 pc) and interactions with companion galaxies. These 
are aspects which have seen significant theoretical progress in recent years. 

We refer to other recent reviews for complementary information on the CGM. 
(2017) provide an excellent overview, but with more emphasis on observational 
properties. (2020) cover the baryon cycle, but also with a focus on obser- 
vations. The recent article by covers both observations and theory, 
but focuses on the more massive halos. Another perspective is provided by [Putman et all 
(2012), who review the observational properties of gaseous halos around the Milky Way and 
other low-redshift spiral galaxies. Readers interested in the physics of the IGM on larger 


scales can refer to the excellent review articles by [Meiksin] (2009) and [McQuinn] (2016). 
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Galactic winds are critical in shaping the CGM but they are a large subject by them- 
selves, so our treatment of winds will be limited in this article. For more on galactic winds, 
readers can refer to [Veilleux et al (2020). Also mostly beyond the scope of this review is 
feedback from active galactic nuclei (AGN; e.g. |Fabian[2012). AGN feedback may have large 
effects on the CGM, but most of the work on AGN feedback has focused on massive (> L.) 
halos and the results typically depend on highly uncertain black hole physics assumptions. 
This is a fascinating topic which would merit much more comprehensive discussion than we 
are able to fit in this article. Among the open questions is whether AGN feedback could 
have important effects on the CGM across a wider range of halos than is often assumed. 


Observations of AGN-driven outflows in dwarf galaxies (e.g., 2019) 
and of the Fermi bubbles in the Milky Way 2010, |Predehl et al.|2020) suggest 


this is a possibility worthy of serious consideration for future research. 


2. COSMOLOGICAL PROCESSES 


Figure [2] sets the stage with an illustration of different physical physical processes found 
in a cosmological simulation following the formation of a main galaxy which by z — 0 will 
have a dark matter halo of mass My œ% 1.4 x 10H Mo (Anglés-Alcázar et al.|2017). It is 
clear that the CGM is a highly dynamic environment intimately tied to the assembly of 
galaxies. Gas flows into dark halos from the IGM that pervades the large-scale structure 
of the Universe. Some of these inflows accrete directly onto the central galaxy, where a 
fraction of the baryons form stars or build up the interstellar medium (ISM). Of the gas 
accreted by galaxies, some or even most is ejected back into the CGM through galactic 
winds before it has time to form stars. These winds can later re-accrete and thus “recycle,” 
potentially a large number of times (more on this in Q.3.3). Moreover, galaxies do not form 
in isolation but rather frequently have satellite galaxies. These satellites can strongly affect 
the CGM of the central galaxy, for example by losing their ISM through tidal stripping, 
ram pressure stripping, or galactic winds of their own (more on this in 92.4). The fractions 
of the total CGM mass originating from fresh accretion, galactic winds, and companion 
galaxies vary as a function of halo mass and redshift, and depend on feedback details, but 
can be comparable to each other especially around ~ L, galaxies (Hafen et al.[2019). Thus, 
all these processes are important to consider in general. 

With the main CGM components identified, the rest of this section reviews in more 
detail our current knowledge of the cosmological processes that shape the CGM. 


2.1. Gas accretion 


2.1.1. Cold vs. hot accretion. Before considering the full complexity of the CGM, it is useful 
to examine the cooling physics of gas in dark matter halos in the idealized approximation of 
spherical symmetry and neglecting feedback processes. Three timescales are important: the 
Hubble time ty = 1/H (where H is the redshift-dependent Hubble parameter), the free-fall 
time in the gravitational potential fg and the cooling time of the gas teoo1. As dark matter 
halos form from gravitational clustering, gas is dragged inward. In the lowest-mass halos 
the gas inflows remain subsonic owing to heating by photoionization by the cosmic ionizing 
background; in those small halos galaxy formation is suppressed (e. g.,[Efstathiou[1992] [Noh] 
& McQuinn|2014). In more massive halos, inflows reach supersonic velocities and are shock- 
heated to a temperature of order the virial temperature, Tyir = (ump/2k)(GMn/Rvir), 
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Figure 2: Different types of circumgalactic gas flows in a cosmological zoom-in simulation 
of a main galaxy with halo mass My z 1.4 x 10 Mo at z = 0. The background gray scale 
shows projected gas density (logarithmically scaled) in a (240 kpc)? volume (physical), while 
the orange stars represent the stellar distribution in different galaxies (the virial radius of 
each dark matter halo is indicated by a dashed circle). The different panels show the system 
at different redshifts. Purple lines indicate streamlines of “fresh accretion,” i.e. gas that 
will accrete onto the main galaxy directly from the IGM by z — 0. Blue lines indicate the 
future trajectories of gas particles ejected from the central galaxy that will accrete back 


^wind recycling" component. Green lines indicate the future trajectories of 


as part of the 
gas particles removed from the ISM of another galaxy that will later smoothly accrete onto 
the central galaxy as part of the “intergalactic transfer" component. Not shown in this 


visualization are galactic winds that escape galaxies and contribute to CGM mass but do 


not re-accrete onto the central galaxy. Adapted from |Anglés-Alcázar et al.| (2017). 


where u is the mean molecular weight (~ 0.6 for an ionized cosmic plasma) and mp is the 
proton mass. For a halo of mass Mn = 10!? Mo at z = 0 (similar to the Milky Way), the 
virial radius Ryir œ 260 kpc and the virial temperature Tue & 6 x 10° K (Barkana & Loeb 


2001). Since Ry "e Mè’ e the virial temperature TV; cx M?! 3. for Z L, halos this gas 


emits in X-rays. The character of gas accretion onto the central galaxy (and of the CGM) 
depends on whether the cooling of the shocked gas is rapid or slow relative to the free-fall 
time. 

Cold accretion. When tcool < tg, the shocked gas rapidly cools and loses its thermal 
pressure support. The cold T ~ 10+ K gas that results tends to fragment and clump, and 
can also form narrow filaments known as “cold flows" or “cold streams" (see Q.14]on cold 
streams and more in $3] about the small-scale properties of cold gas). If unimpeded, e.g. 
by feedback or angular momentum, the cold gas can accrete onto the central galaxy on a 
free-fall time. Since the infall of the cold gas is highly supersonic (relative its internal sound 
speed), a strong shock can form on impact with the central galaxy. 

Hot accretion. When too. > tg, gas cooling becomes a rate-limiting step. Shock- 
heated gas can be supported for an extended period of time ~ teoo1 in the halo potential by 
thermal pressure. In the inner regions, within the “cooling radius" where Lei < tu, there is 
sufficient time for the hot gas to cool and accrete smoothly onto the central galaxy. Absent 
feedback, these cooling regions tend to a steady-state “cooling flow” in which compressional 
heating in the inflowing gas balances radiative losses and Leet 7 £g (e.g., Fabian et al.]1984), 
though in practice feedback processes can modify the flow. 

The different limits corresponding to different regimes of tcooi/t@ are core ingredients 
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of theories of galaxy formation, starting from influential analytic models from the 1970s 
[White & Sec) The implications of 
these limits for galaxy formation as well as the CGM have been the subject of extensive 
investigation since, using analytic and semi-analytic techniques (e.g., 


Somerville & Primack]1999| [Dekel & Birnboim|2006), idealized numerical simulations (e.g., 
Birnboim & Dekel 2003| Fielding et al.12017b | [Stern et al 2020), and detailed cosmological 
simulations (e.g., 2005} |2009b| |Faucher-Giguére et al.12011||van de Voort et al. 


2011 2013). Some ideas are summarized in 82.1.6|though this is still an active 
area of research and (perhaps surprisingly) there is not yet agreement on the effects of cold 


vs. hot accretion for galaxy formation and evolution. 

In the above sketch, we have deliberately been ambiguous about where the cooling 
and free-fall times are evaluated. Modern hydrodynamic simulations as well as observa- 
tions indicate that the CGM can be highly inhomogeneous and consist of multiple phases. 
Therefore, different tcooi/t# limits can be realized in different regions. The physical picture 
is further complicated by outflows from stars and black holes (82.3), as well as additional 
physics such as magnetic fields, thermal conduction, and cosmic rays (83), which imply 


there is in general much more to the CGM than just cooling and gravity. 


2.1.2. Maximum hot gas accretion. To gain further insight into the different modes of gas 
accretion in halos, we consider some analytic results regarding the maximum rate of hot 
gas accretion. Our treatment here follows Stern et al and [Stern et al.] (2020), who 
analyzed the physics of cooling flows in galaxy-scale halos. Although real halos can be much 
more complex and dynamic than idealized cooling flows, this simplified setup allows us to 
develop analytic insights that apply in regions where the gas dynamics is dominated by 
gravity and cooling. These results build on and extend previous on work on cooling flows 
in clusters of galaxies Fabian et al.11984). In clusters it is well 
known that cooling flow models fail to explain the X-ray properties of the ICM. The jury 
is still out as to whether pure cooling flow models can adequately model the CGM of some 
lower-mass systems, since X-ray observations can currently only barely probe the hot gas in 
such halos[!] We do not take a position on this here, but simply use cooling flows as a useful 
baseline solution to gain insight into expected CGM properties before they are modified by 
feedback. 

The setup is a spherically-symmetric dark matter halo in which there is initially a 
pressure-supported, steady flow of gas near the virial temperature. The energy conservation 
equation is v,.d[v2/2 + ye + ®]/dr = —q, where r is the radius, v. is the radial velocity, e 


is the specific thermal energy, y is the adiabatic index of the gas, is the gravitational 
potential, and q is the cooling rate per unit mass. The sum in square brackets is the Bernoulli 
parameter, which is conserved along stream lines in a steady flow. To first approximation 
the first two terms can be neglected for slow inflow and for potentials that are not too far 
from isothermal (such that the specific thermal energy gradient is small), so that d®/dr ~ 
—q/vr = —ngà/ pu,., where A is the cooling function. Since the mass accretion rate M- 
—A4nr? pv,. (the accretion rate is positive when the radial velocity is negative), we have the 


1We note this is plausible since e.g. stellar feedback can in principle act very differently on galaxy 
scales than AGN feedback acts on cluster scales. Moreover, outflows appear to be relatively weak 
around low-redshift ~ L, galaxies such as the Milky Way, so their hot CGM may be reasonably 


well approximated by pure cooling physics (Stern et al./2019). 
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following expression in terms of the gas cooling rate and the radial gradient of the potential: 


ire Anr?n2; A 
do /dr 


At any radius in the halo, there is a maximum steady accretion rate of hot gas, which 
Z tg. At higher 


^2 


is set by the requirement that the density must be low enough that teool 
densities, the rate of compressional heating in the cooling flow cannot balance the radiative 
cooling rate: the gas rapidly cools to «& Tvir. The maximum density can be evaluated using 
te = V2r/ve (where ve is the circular velocity in the potential) and tooo = €/q = pe/nà A: 
` HpUc€ — Mpve 


NH,max(T) We XAr ES XAr £z 0.007 em ^? vioorio Ange, 2. 


where X = 0.75 is the hydrogen mass fraction, vioo = Vc/(100 km s~*), rio = r/(10 kpc), 
and Aan = A/(10 ?? erg cm? s~'). The second equality follows from e z v2, which is 
equivalent to the statement that the gas is at the virial temperature. Using d®/dr = v? Jr, 
the maximum density corresponds to a maximum hot gas accretion rate 


, Anm2vér = E 
Mmax(r) x Fee x3 Mo yr i Viool10À "Zo. 3. 


2.1.3. Virialization of the inner CGM and the threshold halo mass. Equations (2) and 
imply that the maximum hot gas density and accretion rate depend on radius. For gas in 
galaxy halos, the ratio of the cooling time to the free-fall time generally increases from the 
inside out (e.g., [Stern et al.|2020). 'Therefore, the outer parts can be hot and virialized 
(tco /tg > 1) P| while the inner parts cool rapidly and tend toward free fall (tcoo1/t# < 1). 
'The fact that the inner CGM virializes last is important because this defines the time at 
which the boundary conditions of the central galaxy change. 

Cooling flow solutions also reveal an important connection between cooling and whether 
the flow is subsonic or supersonic. In the hot part of the cooling flow, where the temperature 
kTvir ~ myv?2, the sound speed cs ~ V Din ~ Uc. Thus, the free-fall time tẹ ~ r/cs is of 
order a sound crossing time. In this region, the inflow rate is limited by cooling, so we have 
Lei ~ r/|v.|. Combining these results and defining the Mach number M = |v;|/cs, 


in ME 4. 
ig 

In this expression we have omitted a prefactor ~ 1 whose exact value depends on the shape 
of the gravitational potential. It follows from equation that the radius where tooi/t# œ 1 
coincides with the sonic radius Rsonic, where M = 1. The flow is subsonic outside Rsonic 
but supersonic inside that radius. The transition from a subsonic to a supersonic flow has 
important implications for both the physics and observational properties of the CGM. In 
particular, thermal instability is inhibited in the subsonic region of a standard cooling flow 


whereas it can grow faster than the flow time in the supersonic region (Balbus & Soker 


1989). In the supersonic region, large density and pressure fluctuations develop as a result 


By virialized, we mean that a virial-temperature phase is long-lived. Such gas can be sustained 
for longer than either a cooling time or a free-fall time if there is a continuous supply of gas, e.g. 
through accretion from the IGM, because in the cooling flow that develops (if feedback is neglected), 
compressional heating in the accreting gas balances radiative cooling. 
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Figure 3: Critical gas accretion rate for hot gas and the threshold halo mass, indicating 
when halos complete virialization. (a) The critical accretion rate as a function of halo 
mass at z = 0, for different assumed metallicities Z for the gas (Zwzn is the metallicity 
implied by the observed relationship between galaxy mass and ISM metallicity). (b) Similar 
to (a), but with the critical accretion rate converted into a total gas mass by multiplying 
by the Hubble time. The curves all assume that the gas metallicity is consistent with 
the z = 0 mass-metallicity relation, but show the results for different redshifts. In this 
figure, fp = Qbn/OQm is the cosmic baryon fraction and Lon is the fraction of the halo 
baryonic mass that is in CGM gas. (c) The threshold halo mass above which the CGM is 
expected to be completely hot, Mthres, as a function of redshift for different assumptions. 
This corresponds to when the gas mass in the halo is equal to Menta- The fiducial case 
corresponds to a baryon-complete CGM (f = foam = 1) on the mass-metallicity relation, 
and the other curves show how the threshold mass is modified when either the gas mass 
in the CGM or its metallicity is reduced. The threshold halo mass also depends on the 
spin of the gas (via the circularization radius), but this dependence is not shown here for 
simplicity. Adapted from Stern et al. (2020). 


of thermal instability (e.g., Stern et al.|2020), which may have important implications for 


observational signatures as well as how galaxies interact with the CGM via inflows and 
outflows (see (2.1.6). We discuss thermal instability further in where we note that 
if feedback keeps the hot gas close to global hydrostatic and thermal equilibrium, local 
thermal instability can develop if Leed /fe < 10. 

We now address the question of which halos are expected to be virialized. Since a 
given halo can be virialized outside its sonic radius but not inside it, the question is made 
more precise by considering the point at which the CGM becomes entirely virialized, i.e. 
when the sonic radius becomes equal to the radius of the central galaxy. As stressed above, 
whether the CGM can sustain teooi/t# > 1 depends on the complexities of the baryon cycle 
(see Fig. and feedback in particular, as it affects tcoo1 by heating up the gas and by 
changing its density and metallicity. A critical halo mass can however be derived based on 
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simplified assumptions. 

The CGM can be considered to complete virialization when the accretion rate is below 
Mmax all the way to the circularization radius Reire © V2Ryird & 0.05Rvir, where inflowing 
gas becomes supported by angular momentum and which we use as a proxy for the inner 
boundary of the CGM (the spin parameter is defined and discussed further in {2.2}. 
For any given halo, this defines a critical gas accretion rate Merit equal to Minax (oic). 
Approximating the cooling function as A c T~°"Z°®, valid for T ~ 10° — 10’ K and 


metallicities Z 7 0.3Z (Wiersma et al. [2009], [Stern et al) (2020) obtained 
Merit © 0.7 Mo yr! vigoFaoZo.$ 5. 


where v1oo = v./(100 km s 1) is the circular velocity, Rio = Reirc/(10 kpc), and Zo. = 
Z/(0.3Z). The circular velocity and gas metallicity are evaluated at Reirc. The value of 
Merit as a function of halo mass is plotted in the top left panel of Figure |3| for different 
metallicities at z = 0. The bottom left panel of Figure [3] shows Min as a function of 
halo mass for different redshifts, assuming a mass-dependent metallicity consistent with the 
observed mass-metallicity relation for galaxies. 

The critical accretion rate can be translated into a threshold halo mass Mthres by setting 
the total gas mass in the halo Mgas = foam fe My (where fp = Qn/Qm zi 0.16 is the cosmic 
baryon budget and fcaw is the fraction of this budget in CGM gas) to Mesiti: The idea 
is that Meritt is an estimate of the hot gas mass in the halo when virialization completes, 
so the CGM will be fully virialized only for Mi > Minis. The right panel in Figure 
shows Minres as a function of redshift. The solid curve in this panel shows the result for a 
baryon-complete CGM (foam = 1) and other fiducial assumptions. Interestingly, Mthres is 
roughly independent of redshift, staying in the range ~ (1 — 2) x 10!? Mo from z = 0 to 
z = 6. To see why, note that at fixed metallicity, equation implies Merit x v? A pss 
For a matter-dominated universe, at fixed Mn, Reire c Rvir ex 1/(1 + zl, ve ~ Uvir = 
V GMny/ Ry: e (1 + z)i/? (vvir is the virial velocity), and ty œ (1+ gy, Therefore, 
Meritt ex (14 Zoe. which depends weakly on redshift fl 

We note that halos can be virialized substantially below the fiducial threshold mass 
My ~ 10? Mo, for example if the gas metallicity is lower than assumed or if the CGM 
density is below that implied by the cosmic baryon budget. Strong stellar feedback may 
indeed deplete the CGM by large factors in low-mass halos (e.g.|Hafen et al.|2019). In these 
limits, the entire CGM can potentially be virialized in halos of mass as low as ~ 10! Mo, 
or even less. 

The threshold mass derived above based on cooling-flow arguments is similar to thresh- 


old mass previously derived based on the stability of virial shocks (Birnboim & Dekel 2003 
Dekel & Birnboim||2006| see Fig. EI In these derivations, one considers cool gas accret- 


ing supersonically into halos and shocking as the central galaxy is approached in the inner 
regions. The shock is considered stable when the cooling time of the shocked gas is suffi- 
ciently long for its thermal pressure to drive outward expansion of the accretion shock. The 
threshold halo mass derived in this way roughly matches the one derived based on cooling 


?'The bottom left panel of Figure [3] shows that Meritt depends more strongly on redshift for 
low-mass halos. This is because the weak redshift scaling depends on the A œx T? temperature 
scaling, which is only a valid approximation for T ~ 10? — 107 K, where metal lines dominate 
the cooling rate. For lower mass halos, cooling by H and He is important and A has a different 
temperature scaling. 
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Figure 4: Time sequence, from z = 1 to z = 0, showing how the CGM completes virialization 
in a cosmological zoom-in simulation with realistic stellar feedback of a dark matter halo of 
final mass Mi, ~ 10? Mo. Each panel shows gas temperature for a slice through the halo 
center and is labeled by the ratio of the cooling time of the hot, shocked gas to the free-fall 
time evaluated at 1096 of the virial radius. The virialization of the CGM completes when 
tcoo1/te Z 1 in the inner regions, but hot gas is present at larger radii earlier. Before the 
inner CGM virializes, the central galaxy is surrounded by a highly inhomogeneous mixture 
of cold and hot gas. Following virialization of the inner CGM, the halo is filled with hot, 
relatively uniform gas extending to the boundary of the central disk galaxy. Adapted from 
Stern et al. (2021). 


flows because both follow from a comparison of similar cooling and flow timescales. The 
cooling-flow derivation has the advantage of highlighting the fact that the inner parts of 
the CGM "stay hot and virialized" last, which is opposite to the inside-out direction in 
which accretion shocks propagate. The key reason for this difference is that, once hot gas is 
created, whether it will stay hot or rapidly cool in a given region of the CGM is a function 
of the local tcoo1/tg ratio, regardless of the directionality the shock that originally heated 
the gas. On average this ratio increases from the inside out. 

We stress that the cooling flow and virial shock stability treatments are two idealized 
models for gas virialization in halos. The two models provide complementary insights, but 
we do not expect either to perfectly describe the dynamics of the real CGM, which are 
more complex due to time-variable inflows and outflows, as well as strong departures from 
spherical symmetry. In cosmological simulations including realistic feedback, such as the 
FIRE zoom-in simulation of a Milky Way-mass halo shown in F igure [1] it is found that the 
CGM is often first heated out to large radii by shocks due to star formation-driven galactic 
winds, before the theory predicts that pure accretion-driven shocks should be stable[!] As 
a result, the outer parts of low-mass halos can be hot well before cooling times in the inner 
CGM become long enough to sustain a virialized CGM throughout the halo. Although this 
fact has seldom been emphasized in the literature so far, other simulations also find that 
the outer CGM is typically heated to ~ Tyir before the inner CGM is able to virialize (for 
example, this is apparent in temperature profiles of halos from the EAGLE simulations 
analyzed in Correa et al. 2018 and Wijers et al. 2020). 

While it is beyond the scope of this review to discuss detailed observational predictions, 
we note here some possible observational implications of outside-in virialization: (i) The 
cool inner CGM should give rise to a high incidence of strong low-ionization absorbers, such 


^Note that there is evidence that star formation-driven outflows have typical velocities ^ vc (see 
2.3.1), so it can be difficult to distinguish gas that has been shocked-heated by gravitational vs. 
feedback processes, especially after mixing. 
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as Mg II, at small impact parameters. (ii) Since most sight lines to background quasars 
intersect the outer CGM rather than the inner GM (due to area weighting), the presence 
of hot gas at large radii may contribute to the prevalence of multiphase gas inferred in 
observations across a wide range of halo mass. These expectations appear broadly consis- 
tent with the existing quasar absorption line data (e.g., |Tumlinson et al.2017). Outside-in 
virialization would also be good news for observations that aim to detect low surface bright- 
ness rest-UV emission from the CGM, since emission is most sensitive to the inner CGM 
(luminosity L c n?) and these wavelengths probe cool gas (Morrissey et al.|2018). These 
observations can potentially test models of CGM virialization, but the observational signa- 
tures in emission have yet to be worked out quantitatively. 


2.1.4. Cold streams and where we expect them. Because the inflow is subsonic in virialized 
gas, the sound crossing time ts = r/cs is short enough for pressure waves to smooth out 
density fluctuations. As a result, accretion of hot gas tends to proceed quasi-spherically. 
On the other hand, cold gas can clump into much finer-scale structures, as we discuss exten- 
sively in Cosmological simulations predict that cold gas inflows often form filamentary 


structures known cold streams or cold flows (e.g., 
[van de Voort et al.|2011). The panels on the left of Figure |5| shows examples of 
such cold filaments in Mp ~ 10? Mo halos at z = 2 in cosmological zoom-in simulations 
evolved with the moving-mesh code Arepo (Nelson et al.|2016). The simulations shown in 
the figure neglect galactic winds and do not include cooling by metal lines, so these halos 
are significantly above the threshold mass for CGM virialization, which is lower for metal- 
free gas. In this regime, the narrow cold filaments are seen to co-exist in the CGM with 
a volume-filling hot phase. Cold streams have been the subject of much attention because 
in some regimes they could be a primary mode of gas accretion for galaxies (e.g., 
et al.12009a), although whether and when this is the case remains unclear as it depends on 
whether the cold gas survives all the way to the central galaxy during infall through the 
CGM, as well as the efficiency with which hot gas is accreted. 

How can cold streams exist above Minres ~ 107? Mo? pro- 
posed an explanation in terms of the geometry of the large-scale structure, which also 
provides insight into why cold streams in massive halos appear to be a high-redshift phe- 
nomenon (e. g., [Keres et al.[2005). The idea is that halos of different masses are, on average, 
located in different regions of the cosmic web. While low-mass halos tend to be embedded 
in large-scale filaments whose cross sections are larger than halo radii, high-mass halos tend 
to reside at the nodes where large-scale structure filaments meet. Therefore, while the en- 
vironment of low-mass halos is roughly isotropic on the scale of the virial radius, high-mass 
halos are fed by collimated structures. What constitutes a ‘high’ vs. a ‘low’ halo mass 
in this context is determined by the non-linear clustering scale, Mn, i.e. the halo mass 
corresponding to density peaks that become exponentially rare. The key point is that Mii 
increases with time due to the growth of structure, so it is smaller at high redshift. The 
panel on the right of Figure D]shows the halo mass vs. redshift corresponding to lo and 20 
peaks in black dashes. The solid red curve in this panel shows that threshold mass above 
which virial shocks are stable according to b m ao 's analysis (similar to 
the Mthres based on the cooling flow argument outlined in 92.1.3). This plot shows that 
halos of mass Minres are common (< Lol and can be considered low-mass at z S 1 but 
become increasingly rare (> Lol above this redshift. Thus, above z ~ 1 halos more massive 
than Mines are increasingly fed by collimated large-scale structure filaments. The higher 
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Figure 5: (a) Renderings of gas density and temperature in halos of total mass ~ 107? 
Mo at z = 2 in cosmological moving-mesh, zoom-in simulations. These simulations neglect 
galactic winds, so the CGM structure is primarily set by the physics of gas accretion. In 
this regime, a hot virialized CGM fills most of the volume out to beyond the virial radius 
(outermost circles) but dense, clumpy cold streams penetrate deep into the inner halo. 
Adapted from Nelson et al. (2016). (b) Analytic theory for the threshold halo mass below 
which accretion onto galaxies is cold and above which it is hot, assuming a post-shock gas 
metallicity of 0.1Z (solid red). This is similar to the Mtnres mass shown on the right in 
Figure |3| but is derived from a virial shock stability argument rather than cooling flow 
physics. The threshold halo mass is nearly constant with redshift. The inclined purple line 
shows a model for the maximum halo mass Mstream below which cold streams can persist in 
a hot CGM. The black dashed curves show the characteristic mass of newly-forming halos 
vs. redshift, corresponding to 1c and 2c fluctuations in Press-Schechter theory. From Dekel 
& Birnboim (2006). 


densities in filaments, relative to the mean densities in halos, imply shorter cooling times. 
On the other hand, the free-fall times are set primarily by the global mass distribution in 
halos and are mostly unchanged. Their short cooling times enable gas filaments to remain 
cold as they fall into massive halos. The cooling times can be further shortened by com- 
pression of the cold streams by the volume-filling hot phase. The oblique purple curve on 
the right in Figure[5] shows a simple analytic model from [Dekel & Birnboim] (2006) for the 
redshift-dependent maximum halo mass for which cold streams are expected in hot halos, 
Mstream, based on a comparison of timescales taking into account the over-densities of fila- 
ments feeding massive halos. In this model, Mstream © (Mthres/fMni)Mthres, where f ~ 3 
is dimensionless factor calibrated from numerical simulations. 

Although this estimate for the maximum mass of halos expected to contain cold streams 
is a useful guide, it neglects a number of important questions regarding the survival of cold 
gas, especially as it interacts with a hot phase. Whether cold streams survive during infall 
into halos depends on processes, such as shocks and fluid mixing instabilities, that are not 
well resolved in cosmological simulations. We discuss the small-scale physics of cold gas 
survival in much more detail in 

Some early results on cold streams using cosmological simulations were questioned be- 
cause they were obtained using traditional smoothed particle hydrodynamics (SPH) meth- 
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ods, which were shown to suppress fluid mixing instabilities and can lead to the artificial 
survival of cold gas (Agertz et al.[2007] Sijacki et al-[2012). Although the detailed properties 
of cold streams remain uncertain because of the relatively low resolutions in cosmological 
simulations, there is currently a broad consensus between different modeling methodolo- 
gies that the existence of cold streams is a robust theoretical prediction. Cold streams are 
found not only in cosmological simulations evolved with modern SPH codes, which have 
been improved to more accurately capture mixing, but also in simulations using adaptive 
mesh refinement (AMR), moving mesh, and mesh-free codes (for a comparison including 
several of these methods, see [Stewart et al.[2017). It is also noteworthy that cold streams 
are found in simulations that vary by orders of magnitude in resolution, ranging from large 
cosmological boxes to cosmological ‘zoom-in’ simulations focusing on individual halos (e.g., 
[2016]. Nevertheless, it is important to keep in mind that cosmological 
simulations still fall short of capturing all the physics relevant to cold gas formation and 
survival, so the theory of cold streams could still evolve substantially. Approaches that 
incorporate insights from small-scale studies will play an important role going forward (see 
84). 

On large scales, interactions with galactic winds and with satellite galaxies can also 
modify the properties of cold streams. For example, galactic winds (including winds blown 
by dwarf galaxies embedded in cold streams) can *puff up" the cold gas distribution 
2015). The increased cold gas cross section in halos due 
to winds and galaxy interactions (see also has important implications for observables, 
such as the predicted cross section for Lyman limit absorption. 


2.1.5. Absorption and Lyo emission from cold streams. Cold streams are of interest as 
observables in the CGM owing to their relatively high densities and their temperatures 
T ~ 10* K. In absorption, cold streams are predicted to manifest as HI absorbers with 


columns in the range Nm ~ 10'° — 107° cm~?, corresponding to Lyman limit systems 


(LLSs) and partial LLSs (e.g., |Fumagalli et al.12011b||2014| |Faucher-Giguére & Keres12011 
Faucher-Giguére et al.|/2015 2017). Cold streams may in fact dominate the 


incidence of these strong absorbers at most redshifts where they are observed (e.g., [van de] 
2012) and metal-poor LLSs have been interpreted as detections of cold streams 
infalling from the IGM which have not yet been significantly enriched by feedback processes 


(e.g., Ribaudo et al.]2011| |Fumagalli et al.12011a). 


In emission, cold streams may be important in explaining spatially extended structures 
known as Lya halos (e.g., |Steidel et al.[2011) or the more extreme Lya blobs (e.g., 
et al. Rescue Un Cro RE Ba One possibility is that gravitational 
energy is released as Lya “cooling radiation” during the infall of cold streams (c.g., [Dijkstra] 
[& Loeb[2009). A simple estimate shows that cooling radiation could in principle be very im- 
portant. Let Meas be the gas accretion rate in the halo and A the difference in gravitational 
potential as gas falls from the IGM down to the inner halo. Assuming an NF W potential 
with concentration c = 5 and a gas accretion rate Meas = fo Miot, 
where Mo, is an average total mass accretion rate following [Neistein & Dekel (2008), the 
cooling luminosity L”! ~ fa ef Mgas|A®| ~ 4 x 10% erg s7! fa ea Mis®[(1+2)/4]*°, where 
faes is an efficiency factor quantifying how much of the gravitational energy is released in 
the Lya line and Mi; = M,,/(10'* Mo) (see the appendix in [Faucher-Giguere et al.[2010). 
This luminosity is comparable to observed Lya halos. 

However, the temperature of cold streams puts them on the exponential part of the Lya 
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emissivity function. Namely, the Lya emissivity powered by collisions is e£?! = Ca(T)nuine, 
where C is the collisional excitation coefficient, nu; is the neutral hydrogen number den- 
sity, and ne is the free electron number density. The collisional excitation coefficient scales 
as Ca x T-/? exp (—T. /T), where Ty = hva/k © 1.2 x 10° K and va is the Lya frequency. 
This exponential dependence on temperature makes theoretical predictions for cooling ra- 
diation highly uncertain (Faucher-Giguére et al. 
et al.[2020b). In simulations, the predictions are sensitive to the numerical methods used to 
model the hydrodynamics (because of the importance of fluid mixing instabilities and weak 
shocks) and radiation (because it alters the ionization structure and photoionization also 
heats the gas). The structure of turbulent mixing layers at the boundaries between cold 
and hot gas, discussed in is relevant as it may be where much of the energy dissipation 
occurs, but these layers are not resolved in cosmological simulations. 

Alternatively, extended Lya emission can be powered by recombinations following ion- 
ization by stars or AGN. These recombinations can occur either in the ISM (HII regions) 
or, for ionization radiation that escapes galaxies, in the CGM. In the case of Lyo photons 
produced within galaxies, diffuse halos can be formed by resonant scattering with neutral 
hydrogen in the CGM (e-g. [Dijkstra et aT [2006] [Gronke et aL [2015]. For reference, the Lya 
emission powered by stellar radiation in HII regions L” ~ 104% erg st fo,escSF Rio, where 
fa,esc is the fraction of Lya photons that avoid destruction by dust and escape the medium 
and SFRio = SFR/(10 Mo yr^!) (e.g., [Leitherer et al.|[1999). This is comparable to the 
Lya luminosity of cooling radiation, which in part explains why it has been difficult to 
unambiguously identify what powers observed sources (scattering can in principle be tested 
using polarization; [Dijkstra & Loeb|2008). In the case of ionizing radiation that escapes 
galaxies, Lya photons can be produced in the CGM via fluorescence, i.e. recombination 
emission powered by ionizing photons absorbed in the halo 


[et al |[2010). The Lya emissivity from recombinations e" = a. ,recani(T)nunne, where 
fa,rec is the average number of Lya photons produced per recombination (fa,rec % 0.68), 


anı(T) e T "HI is the recombination coefficient, and num is the ionized hydrogen number 
density. While recombinations are not as sensitive to temperature as collisional excitation, 
the recombination emissivity is sensitive to gas clumping (the emissivity is proportional to 
the clumping factor C — (n?)/(n)?). This dependence on the clumping factor has been 
used to infer unexpected small-scale structure in the cold gas in the halos of some luminous 
Lya blobs (Cantalupo et a1.[2014] [Hennawi et al.J2015). This has led to the proposal that 
the CGM could be filled with a “fog” or “mist” of tiny but high-density cold clouds; the 
physics of these tiny clouds is covered in 

Galactic winds can also power extended emission by depositing mechanical energy into 
the CGM, which can then be radiated away (Taniguchi & Shioya[2000] [Sravan et al.[2016). 
Even if the ultimate energy source for extended emission (whether it be radiation or me- 
chanical energy from stars and/or AGN) originates from galaxies, cold streams may be 
important to explain Lya emission on halo scales. This is especially the case in massive 
halos exceeding the threshold mass Mthres ~ 10!” Mo, above which the volume-filling phase 
is expected to be hot. If all the halo gas were hot, most of the emission would be expected 
to come out in X-rays. Cold streams and other cold gas structures in halos, such as a 
possible cold ‘fog’ (§3.3), can scatter Lya photons that escape galaxies or otherwise ensure 
that a significant fraction of the energy deposited into the CGM is radiated in Lya rather 
than in higher energy bands. 
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2.1.6. Effects of CGM virialization and accretion mode on galaxies. Much of the interest in 
the CGM is rooted in the presumption that the physics of gaseous halos plays an important 
role in the formation of galaxies. In particular, there is broad but indirect observational 
evidence that CGM virialization is important for galaxy evolution. The characteristic lu- 
minosity of galaxies, L, (above which the galaxy stellar mass function is exponentially 
suppressed), corresponds to a roughly constant halo mass My, ~ 10!? Mo, only weakly 
dependent on redshift. This is also the halo mass scale above which the fraction of galaxies 
that are quiescent rises above ~ 50% (e.g., [Behroozi et al.|2019). In the last few years, 
observations of spatially resolved galaxy kinematics have suggested that the L, mass scale 
is consistent with the emergence of large disk galaxies (e.g.,|Tiley et al.2021). This mass 
scale, termed the “golden mass” by some authors (e.g., [Dekel et al.[2019), is similar to the 
halo mass at which the CGM is theoretically expected to complete virialization (see Fig. 
3). 

Despite the substantial evidence that CGM virialization correlates with major changes 
in galaxy properties, whether and how CGM physics affect galaxy evolution remains an 
active area of research, with basic questions still the subject of debate. We summarize 
below some ideas that have been proposed for how CGM processes could affect galaxy 
evolution for L ~ L, galaxies, and which in our view deserve deeper investigation: 

A quasi-isotropic, hot CGM is necessary for effective preventative feedback. 
There is a broad consensus that in order to explain the observed population of “red and 
dead” galaxies at the massive end, it is not sufficient for feedback to eject gas from galaxies. 
There must also be “preventative feedback” which prevents halo gas from cooling and 
raining onto galaxies at overly high rates (e.g., [Croton et al.][2006). 
In the most massive halos, this feedback is often assumed to come from jets powered by 
AGN, but wider-angle winds powered by either AGN or supernovae can play a role (Type 
Ia supernovae can be energetically important in ellipticals with old stellar populations; e.g. 
[Voit et al.[2015]. An idea often discussed in this context is that preventative feedback only 
becomes important after most of the CGM has become hot and quasi-isotropic (e-g., [Keres] 
et al.[2009a). This is because only in this limit can feedback keep the gas hot. In contrast, 
when there are massive inflows of clumpy or filamentary cool gas, the smaller geometric 
cross section of the inflows strongly reduces the efficiency with which feedback couples to 
accreting gas. 

Pressure fluctuations change at the order-of-magnitude level at inner CGM 
virialization. Whether the CGM is virialized or not also changes the boundary conditions 
of the central galaxy. Both idealized simulations and cosmological 
simulations show that when the inner CGM virializes, there is a change 
from order-of-magnitude thermal pressure fluctuations in the gas around the galaxy (prior to 
virialization) to a roughly uniform pressure (after virialization; see Fig. D. Large pressure 
fluctuations in the inner CGM create paths of least resistance through which feedback can 
more easily expel gas from the galaxy. Thus, we may expect that large-scale galactic winds 
will be stronger and reach farther out before the CGM virializes. There is some evidence 
from galaxy formation simulations with resolved ISM physics that star formation-driven 
outflows are suppressed when the inner CGM is virialized, such as around Milky Way-like 
galaxies at z ~ 0 (e.g. Muratov et al.[2015] [Stern et al [20212]. Large pressure fluctuations 
in the inner CGM may also make it difficult for the ISM to reach a statistical steady state, 
which could result in highly time-variable (or *bursty") star formation rates (e.g., 


fetal. 922). 
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CGM virialization changes the buoyancy of supernova-driven outflows. 
and [Bower et all proposed a related but different effect of CGM viri- 
alization on outflows. These authors suggested that supernova-inflated superbubbles are 
buoyant in the CGM prior to virialization, so that outflows can be “lifted” in the CGM by 
buoyancy forces, but that the bubbles would cease being buoyant once a hot CGM develops. 
These authors argued that stellar feedback would therefore become ineffective at expelling 
gas once the CGM virializes. They furthermore hypothesized that this would lead to the ac- 
cumulation of gas in galaxy centers, which would allow nuclear black holes to start growing 
more rapidly. If correct, this mechanism would represent another connection between CGM 
virialization and AGN feedback. Similar phenomenology regarding accelerated black hole 
feeding starting around L,, found also in other simulations, has however been attributed 
by other authors to changes in star formation-driven outflows either due to confinement 


by gravity or to the pressure fluctuations effect mentioned above (e.g., 
Byrne et al. [2022]. 

Hot accretion promotes the formation of thin disks by making angular mo- 
mentum (AM) coherent. Recently, reported evidence from cosmo- 
logical simulations that hot-mode accretion promotes the formation of galaxies with thin 
disks, such as observed in low-redshift Milky Way-like galaxies. The basic idea is that gas 
from large-scale structure enters dark matter halos with a broad distribution of specific 
AM (sAM). When the gas falls in toward the galaxy as cold clumps or filaments, spatially 
separated gas parcels are causally disconnected. In this regime, the cold gas reaches the 
halo center supersonically with a still-broad sAM distribution and tends to form stars in 
irregular or thick disk morphologies. On the other hand, when the gas accretes onto the 
central galaxy in a smooth, subsonic cooling flow, the sAM distribution becomes coherent 
(i.e., narrow) before accretion onto the galaxy, and stars form in a thin disk configuration] 

The role of the gas accretion mode in determining the morphology of galaxies is an 
example of how there is not yet a consensus on the role of CGM physics in galaxy formation. 
While the recent work mentioned in the previous paragraph highlights the role of hot 
mode accretion in the formation of thin disks, a substantial body of work has instead 
emphasized the role of cold streams in feeding massive disks at high redshift (e.g., 
et al.|2009b). These results are not necessarily inconsistent because the disks in the massive, 
high-redshift regime are highly turbulent and geometrically thick. More work on the role 
of the gas accretion mode on the formation of disk galaxies will be important, including 
special attention to how results vary as a function of halo mass and redshift. 

Despite the plausible causal CGM mechanisms summarized above, we must stress it has 
proved challenging to disentangle whether CGM changes cause changes in galaxy properties, 
or whether changes in CGM and galaxy properties simply correlate. For example, analytic 
arguments suggest that the mass scale of CGM virialization is similar to the mass scale 


where supernova-driven outflows become confined by gravity (e.g., |Lapiner et al.||2021 
Byrne et al.|/2022), so it is possible that outflows are suppressed around the same time 


as the CGM virializes but that neither change drives the other. We conclude that more 
research is needed to firmly establish the roles the CGM plays in galaxy evolution. 


Dit is likely that the net result depends not only on how the gas accretes but also on feedback, 
because absent feedback we would expect a thin gas disk to eventually form as a result of dissipa- 
tion, even if the sAM distribution is not initially coherent (as in other astrophysical settings, e.g. 
protoplanetary disks that form in turbulent molecular clouds). 
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2.2. Angular momentum 


We now expand on what is known about the AM content and exchange processes in the 
CGM. The motivation for this is twofold. First, in the standard cosmological picture, 


galaxies inherit their AM from gas accreted from their host halos (e.g.,|Fall & Efstathiou 
1980 1997 1998). In this picture, the AM of halos is first 


acquired via gravitational torques during structure formation (e.g., 1969 
1984), although the AM of a particular halo fluctuates substantially over time due to mergers 


(e.g., [Vitvitska et al.[2002]. While on sufficiently large scales the baryons are expected to 
have AM properties similar to the dark matter, hydrodynamic forces and feedback processes 
experienced by the baryons during galaxy formation can potentially strongly affect both 
the AM content of the CGM as well as of galaxies. Second, observations indicate that in 
many systems the CGM has substantial rotation (e.g., 
‘Hodges-Kluck et al.12016) and we would like to understand these CGM observations. 

A useful basis to understand the AM properties of the CGM is to start with scalings 
derived from dark matter-only simulations. Dark matter halos can be characterized by a 
dimensionless spin parameter A = Leier: where J is the total AM inside a sphere of radius 


R containing mass M, and V = ,/GM/R is the circular velocity at R (Bullock et al.|2001). 
We adopt the standard choice of setting R to the virial radius. With these definitions, 


numerical simulations have shown that dark matter halos have a lognormal distribution of 
spin parameters, nearly independent of mass and redshift, with a median A z 0.035 (e.g., 
Bett et al. TOU SUED Paina i eu UE um E Serieel. Defining the sAM j = J/M and noting that 

Ryir x My (for halos defined to have constant over-density relative to the mean 
matter ee ped j x mM?’ 3 / VÀ +z. Thus, at fixed redshift j increases with halo mass 
ex M? ?; while at fixed halo mass j increases with time c 1/yT F z as redshift decreases. 
While these scalings apply to dark matter-only simulations, (2020) find 
that the same scalings roughly describe the CGM AM trends with halo mass and redshift in 
the IllustrisTNG hydrodynamic simulation, which includes feedback from galaxy formation. 
Assuming that the sAM of the CGM is comparable to that of the dark matter halo (though 
with significant differences discussed below), the small spin parameters A ~ 0.035 imply that 
AM support is negligible in most of the CGM, becoming only important at a circularization 
radius Reire © V2ARvir ~ 0.05 Ry. The small contribution of rotation to the support of 
halo gas has been confirmed by a systematic analysis of different support terms in EAGLE 
simulations (Oppenheimer}2018). 

Next, we summarize some key results concerning how the CGM AM relates to other 
components, including the dark matter halo and the central galaxy. We also review physical 
explanations for the differences found between the different components. We refer to Figure 
[6] for some quantitative results on spin parameters of the gas and the dark matter around 
simulated galaxies with a halo mass My ~ 10!? Mo at z ~ 2. 


2.2.1. AM in the CGM vs. the dark matter halo. Within the virial radius, CGM gas 
has systematically higher sAM than the dark matter. Interestingly, this is the case even 
in non-radiative simulations, so part of the difference can be attributed to hydrodynamic 
interactions that do not involve cooling (e.g., [Zjupa & Springel[2017). For example, when 
two halos merge, ram pressure will cause the gas mass to become offset from the dark 
matter. Since the simulations also predict that the gas and dark matter spins are on 
average misaligned by ~ 35°, mergers could on average spin up gas more than the dark 
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Figure 6: (a) Radial profiles for the spin parameter of different mass components in 29 
simulated halos analyzed in Danovich et al. (2015). The halo mass ranges from 10/5? 
to 10/7? Mo at z ~ 2. The total gas mass is divided between cold and hot components, 
corresponding to T < 10° K and virial-temperature gas, respectively. The symbols show the 
mean spin in each shell, while the error bars show the standard deviation. The dashed line 
shows a slope of unity. Several trends are apparent: (i) for all components, the mean spin 
parameter increases from the inside out; (ii) the spin parameter of the gas is systematically 
higher than that of the dark matter; and (iii) the higher spin in the gas is primarily driven by 
the cold gas. (b) Evolution of the mean spin parameters inside the virial radius as a function 
of scale factor a. Here the cyan “cold gas+stars” component corresponds to < 0.1Rvir, a 
proxy for the evolution of galactic baryons. Averaged within the virial radius, the sAM of 
the cold gas exceeds that of the dark matter by a factor ~ 2 — 3. (c) Comparison of the 
radial profiles of torques on cold halo gas due to gravity vs. torques due to gas pressure, 
for simulation snapshots ranging from z ~ 1.6 to z ~ 3.2. The q parameter is proportional 
to |r|/|l|, showing that gravitational torques are more important for the cold gas. Adapted 
from Danovich et al. (2015). 


matter. 

When the total CGM mass is divided between cold (T « 10? K) and hot (virial- 
temperature) gas, it becomes clear that higher gas sAM relative to the dark matter is 
primarily driven by the higher sAM of the cold gas, which can exceed that of the dark 
matter by a factor ~ 2 — 3 within Ryir. This indicates an important role for gas cooling. 
analyzed the torques experienced by the dark matter and cold gas 
as they approach the virial radius of halos, and argued that the excess quadrupole moment 
of the cold gas relative to the dark matter could explain the additional sAM acquired by 
the infalling cold gas as a result of more efficient tidal torquing. Specifically, the elongated, 
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thin-stream geometry of the infalling cold gas (relative to the thicker dark matter distribu- 
tion) enables the cold gas to acquire AM more efficiently via tidal torques. In other words, 
cold streams are where cold gas gets extra torque. 

An additional timescale effect contributes to the higher sAM of the cold gas relative to 
the hot gas. Whereas the cold gas can accrete onto the central galaxy on a free-fall time, 
the hot gas is supported in the halo by thermal pressure for at least a cooling time. Thus, 
while much of the cold CGM has typically only recently entered the halo, the hot CGM 
has been built up over a longer period in the past. Since the sAM of matter accreting from 
large-scale structure increases with time, this timescale effect alone tends to enhance the 
sAM of the cold gas relative to the hot gas. The increasing sAM of matter accreting from 
large scales likely explains, at least in part, why the spin parameter of the different halo 
components increases systematically with radius in Figure [6] 

Galaxy formation feedback can also increase the sAM of the CGM relative to the dark 
matter. Namely, the ejection of gas from galaxies by star formation or AGN-driven outflows 
occurs primarily from the inner parts where the baryons have relatively low sAM 
[Springel[2017). If sufficiently strong, feedback can eject some of the low sAM gas not only 
from galaxies but from halos altogether. The preferential ejection from halos of low sAM 
gas also tends to enhance the sAM of the remaining CGM relative to the dark matter. 

Although the detailed quantitative predictions depend on the simulation code, including 
the feedback model, the high sAM of the CGM relative to the dark matter (especially for 
the cold gas) appears robust, as similar results have been found in cosmological simulations 
using different codes (e.g., Stewart et al. 2017; DeFelippis 2020) and for halos in different 
mass ranges (e.g., Oppenheimer 2018). The high sAM of cold gas can produce extended 
rotating structures that have sometimes been called “cold flow disks" which may have ob- 
servational signatures in low-ionization absorption systems co-rotating with central galaxies 


(e.s. Stewart et aL. 2017] 2013). 


2.2.2. AM in the CGM vs. the central galaxy. The relationship between the sAM of galaxies 
and that of their host halos merits some comments. On the one hand, observational studies 
(e.g., and many numerical simulations (e.g., 
2022) find that on average the size of galaxies scales with the virial 
radius of the dark matter halo, with a normalization roughly consistent with that expected 
if the sAM of the galaxy is comparable to that of the dark matter halo. Moreover, it is 
found in some simulations that at fixed stellar mass, halos with larger spin parameters 
on average host larger galaxies (Rodriguez-Gomez et al.|2022). On the other hand, some 
simulations that reproduce the average trend between galaxy size and halo size indicate 
that on a halo-by-halo basis, the spin parameter of the central galaxy is barely correlated 
with the spin parameter of the dark matter halo, when these are measured at the same final 


time (Garrison-Kimmel et al.|2018 2019). It is also noteworthy that the AM 


vector of the CGM is in general misaligned with that of the stars in the central galaxy by 
large angles ~ 30 — 60° (DeFelippis et al. [2020). 

These results could be understood if, to first order, the sAM of galaxies scales with sAM 
of the host dark matter halo when the galaxy is assembled but there is order-unity scatter 
introduced between the sAM of the baryons and of the dark matter over time. For example, 
showed that the spin parameter of a dark matter halo fluctuates by 
factors up to ~ 2—3 due to halo mergers. Because of the very different spatial distributions 
of matter, we expect galaxies to be torqued differently during mergers compared to the 
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much larger halos. The partial decoupling of the spin parameter of galaxies from their 
host halos over time is consistent with the finding of [Garrison-Kimmel et all that 
the stellar morphology and kinematics of simulated Milky Way-mass galaxies are poorly 
correlated with the properties of the final dark matter halos (including spin), but that the 
galaxy properties correlate much better with dark matter halo properties evaluated at the 
time when 50% of the stars had formed. 

Processes other than mergers can also contribute to differences between the AM content 
of galaxies and their halos. One likely relevant factor is that the minority of baryons that 
end up in galaxies, relative to the cosmic baryon fraction, is not necessarily representative 
of the majority of halo baryons. The fraction of baryons found in galaxies peaks at ~ 0.2 
for Milky Way-mass halos and is as low as < 10^? for dwarf galaxies and for central galaxies 
in massive clusters (e.g., [Behroozi et al.[2019]. Another possibility is that the sAM of gas 
accreting through the CGM is not strictly conserved but rather experiences exchanges with 
other components. 


2.2.3. Gravitational vs. gas pressure torques. We do not yet have a detailed understanding 
of AM transport in the CGM, but several mechanisms can contribute. |Danovich et al. 
(2015) decomposed the total Lagrangian torque on gas elements into three components: 


T = dl/dt = T + TP +Ts, where l is the AM vector, To = —pr x V9 is the torque 
due to gravitational forces, Tp = —r x VP is the torque due to pressure gradients, and 
Ts = —lV -v corresponds to viscous stresses. The viscous stress term is negligible in the 


ideal hydrodynamics limit. The bottom right panel in Figure[6|compares log q radial profiles 
for gravitational and pressure torques acting on cold gas for simulated halos at z ~ 1.6 — 3.2 
from [Danovich et al. (2015), where ger « |rs,»|/|l|. The results indicate that the torques 
on cold streams are dominated by gravity rather than gas pressure. These gravitational 
torques are sourced by anisotropies in the matter distribution, ranging from large-scale 
structure to central disks, which tend to align the infalling gas. It would valuable to extend 
this kind of analysis to other regimes in the future. For example, the relative importance of 
gravitational torques vs. pressure torques could be very different for hot gas, which tends to 
be more spherical in geometry and in approximate hydrostatic equilibrium throughout the 
halo. It would also be worthwhile to quantify the effects of magnetic fields on AM transport 
in the CGM. Magnetic fields play a key role in transporting AM in accretion disks around 
young stars and black holes, but their effects on AM exchanges in the CGM has not yet 
received attention to our knowledge. 

After baryons are accreted by galaxies, gravitational torques due to asymmetric features 
in the potential (e.g., spiral arms, bars, or massive clumps) can strongly affect the AM 
distribution within galaxies, such as by forming central bulges (e.g.,|Shlosman et al.]1989). 
'This would also contribute to differences between the AM of galactic components relative 
to what may be expected from strict conservation of AM inherited from the halo. 

Overall, AM acquisition and exchange processes in the CGM remain relatively under- 
studied and more work on this topic would be highly valuable. 


2.3. Galactic winds 


Galactic winds are commonly observed and are an essential ingredient of modern galaxy 
formation theories. In most galaxies, these outflows are understood to be primarily driven 
by energy and/or momentum produced by massive stars, including via supernovae (e.g., 
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Chevalier & Clegg|1985) and/or radiation pressure (e.g., [Murray et al.|2005). In galaxies 


with luminous active galactic nuclei (AGN), galactic winds can also be powered by accretion 
onto massive black holes (e.g., and references therein). 
Here we focus on outflows powered by star formation. In current models, these outflows are 
critical to suppress star formation in galaxies up to ^ L,, either by ejecting gas from the ISM 
before it has time to turn into stars or by preventing CGM gas from accreting onto galaxies 
in the first place [Naab & Ostriker][2017). However, because 
they originate on small scales and their driving mechanisms are not yet fully understood, 
the properties of galactic winds remain highly model dependent. Therefore, we limit our 
discussion below to general concepts and results which are useful to understand the impact 
of galactic winds on the CGM, and vice versa, rather than the detailed predictions of specific 
models. 

F igure [7| summarizes some key properties of galactic winds in simulations. One salient 
feature is that galactic winds are multiphase. This multiphase structure is clearly ob- 
served in the protypical example of the galactic wind driven by the M82 starburst galaxy 


(Strickland & Heckman|2009) and is also predicted by several models. In the models, the 


multiphase structure typically consists of a hot fluid heated by supernova explosions with 


a spectrum of embedded cool clouds (e.g., Schneider et al. 


& Bryan|2022). However, winds accelerated by gentler processes such as radiation pressure 
or cosmic rays could be cooler overall (e.g., [Booth et al.|2013); the 
dominant driving mechanisms and wind properties could well vary with galaxy mass and 
redshift. Much of our discussion of the physics of multiphase gas in including processes 
that govern cold cloud growth and survival, is relevant to galactic winds. We discuss a global 
thermal instability in winds in 93.1.3] and we discuss cloud-wind interactions extensively in 


2.3.1. Bulk scalings. Although the gas in galactic winds exhibits a range of velocities, den- 
sities, and temperatures (even in individual galaxies), there is some evidence that the mean 
(or median) velocity vy scales linearly with the circular velocity of the galaxy (vc). This ve- 
locity scaling is predicted, for example, in FIRE simulations in which galactic winds emerge 
from the energy and momentum injected by multiple stellar feedback processes (including, 
SNe, stellar winds, and radiation pressure) on the scale of individual star-forming regions 
(Muratov et al.][2015]. In large-volume simulations in which the generation of galactic 
winds is not resolved but the wind properties are instead prescribed, it is also found that a 
wind injection velocity proportional to ve can produce a reasonably good fit to the observed 
galaxy stellar mass function (e.g.,|Davé et al. 2014). While we do 
not understand this scaling in detail, we can heuristically reason why it may emerge from 
the self-regulation of stellar feedback (e.g., [Murray et al.|[2005]. Namely, v. scales with 
the escape velocity in the potential, so much slower outflows would quickly fall back onto 
galaxies, strongly suppressing their net effect. On the other hand, much faster outflows 
would easily escape halos and halt galaxy formation. 

The scaling with circular velocity can be used to derive how the mass outflow rate My, 
scales in different limits. When the wind is energy-driven, the product My v2, is fixed, so 
the mass outflow rate scales as My œ 1 /v2, ox 1/v2. Similarly, when the wind is momentum- 
driven, the fixed product is My vw, so My œx l/vw c 1/v.. Since the feedback energy scales 
with the star formation rate, these scalings are often expressed in terms of the mass loading 
factor nu = My /SFR. An example of an energy-driven wind is a hot, supernova-driven 
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Figure 7: (a) Simulation of a galactic wind driven from a disk galaxy by supernovae. The 
density rendering highlights the multiphase nature of the outflow, with a low-density fluid 
heated by supernova explosions escaping the galaxy with a highly structured spectrum 
of cool clouds embedded in it. From Schneider et al. (2020). GAAS. Reproduced with 
permission. (b) Comparison of wind mass loading factors (defined as rw in the main text) as 
a function of dark matter halo mass in different cosmological simulations and semi-analytic 
models (labeled in the legend) at z = 0. The solid curves correspond to outflow rates from 
the galaxy, while the dashed curves correspond to outflow rates from the halo (at one virial 
radius). The outflow rate from the halo can be larger than from the ISM because of gas 
entrainment in the CGM. From Mitchell et al. (2020b). (c) Fraction of the final stellar mass 
contributed by different accretion channels as a function of halo mass at z — 0, for different 
cosmological simulations. The accretion channels are defined similarly as in Figure [| and 
correspond to gas that formed stars after directly accreting onto a galaxy from the IGM 
(“first infall”), gas that recycled in winds before turning into stars, and gas that transferred 
from one galaxy to another before forming stars. From Mitchell et al. (20202). 


outflow in which radiative losses are negligible. An example of a momentum-driven wind 
would be one driven by radiation pressure, in which the momentum of photons is transferred 
to the gas but thermal energy plays a negligible role in the outflow expansion. The top 
right panel in Figure[7] compares mass loading factors measured from different cosmological 
simulations and semi-analytic models as a function of halo mass. 

We note that, because mass loading can occur both in the ISM and in the CGM (see 
below), while the energy and/or momentum injection is concentrated in the galaxy, the 
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energy and momentum loading factors gg = Bw / Eteedback and np = pw/DPteedback are of- 
ten more robust model predictions. Here, the subscript ‘feedback’ refers to the energy or 
momentum injected in the ISM by feedback processes while the subscript ‘w’ refers to the 
energy or momentum escaping in a wind. The energy loading factor gg can be < 1, e.g. 
when the majority of the energy from SNe is radiated away in the ISM before wind break 
out (e.g., Fielding et al.]2017a). 


2.3.2. Entrainment of CGM gas by galactic winds. The properties of outflows can change 
greatly as they expand into the CGM. As outflows expand, they are decelerated by gravity 
as well as by entrainment of CGM mass! The entrainment of CGM gas modifies the mass 
outflow rate, as well as its chemical composition by mixing gas recently ejected from the 
galaxy with ambient CGM. CGM entrainment can be very important: for example, [Muratov] 
and [Mitchell et al.] showed that mass outflow rates at the virial radius 
can be dominated by entrained gas, in the FIRE and EAGLE simulations, respectively (for 
EAGLE, this is shown by ‘halo scale’ mass loading factors that are larger than ‘galaxy 
scale’ loading factors in Fig. H In other simulations, such as IllustrisTNG, the entrained 
CGM mass is less important relative to the gas directly ejected from the ISM 
2019), again underscoring the model dependence of outflow results. Since the metallicity of 
the CGM is generally lower than that of the ISM, entrainment tends to dilute the outflow 
metallicity. These effects imply that it is critical to specify where outflow properties are 
measured (such as at what radius) when comparing model predictions to observations, or 
different models to one another. 

As they sweep up CGM, galactic winds can affect the properties of gas accretion in halos. 
For example, outflows may push out infalling gas and prevent some of it from accreting onto 
galaxies (e.g., [Tollet et al.[2019). Outflows may also affect the survival 
of cold streams, drive turbulence or inject heat into halo gas. Thus, it should be borne in 
mind that galactic winds are likely to modify some aspects of our simplified discussion of 
gas accretion in halos (82.1) in model-dependent ways. 


2.3.3. Wind recycling. An important property of galactic winds is that some or most of 
their mass can recycle, i.e. re-accrete onto galaxies (see Fig. D This implies that, in an 
instantaneous sense, some CGM gas that is observed to be infalling onto galaxies may have 
been previously part of a wind (e. g, [Hafen et al.[2020). A phase change may occur as winds 
recycle, e.g. if a hot wind cools and cold clouds rain back onto the galaxy, but this does 
not necessarily occur if gas is ejected cold from the galaxy, as in some momentum-driven 
wind models. Recycling has also been shown to be very important in an integrated sense in 
shaping the galaxy stellar mass function, as was shown for example in the pioneering study 


of wind recycling by [Oppenheimer et al (2010). The fraction of wind mass that recycles in 
depends on galaxy mass, redshift, and on the feedback model 
can be more than half and up to ~ 1 in some simulations (e.g., 
Auglés-Alcázar ct al 2017). 

Useful concepts to characterize wind recycling include the distribution of recycling times 


and the distribution of the number of recyclings. Long recycling times mean that, after being 
ejected in a wind, a gas element spends a long time outside galaxies before being re-accreted. 


Din {3.2| we will discuss the entrainnment of cold clouds in hot winds. Here, the entrained CGM 
mass can be volume-filling hot gas as well as cold gas. 
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A given gas element can in general be recycled many times. In the FIRE simulations, the 
star formation histories of dwarf galaxies are highly time-variable, and gas elements can be 
ejected then re-accreted up to ~ 10 times by redshift zero. The multiple cycles of wind 
ejection and re-accretion in these dwarf galaxies may be an important factor driving the 
burstiness of star formation predicted by the simulations (Anglés-Alcázar et al.[2017). 

As we discuss below in another form of wind recycling occurs when gas ejected by 
one galaxy re-accretes onto another galaxy. 


2.4. Satellite galaxies 


In this review, we focus primarily on physical processes operating in the CGM of central 
galaxies, i.e. main galaxies at the center of dark matter halos. The CGM of satellite galaxies 
can be affected by additional effects and does not separate cleanly from the CGM of the 
central galaxy they orbit. In this section, we briefly list some of the ways in which satellite 
galaxies can affect the CGM of the central galaxy. 

Most directly, gas that remains bound to satellite galaxies (e.g., satellite ISM) can give 
rise to strong absorption features in the spectra of background sources. If the satellite 
is faint, the satellite may not be detected in emission and the absorption features can 
be mistakenly attributed to the CGM of the central galaxy. Similarly, satellites could 
contribute to the spatially extended emission that sensitive experiments aim to detected 
from the CGM. 

Gas originally belonging to satellites can also be lost and incorporated into the CGM 
of a central galaxy by several different processes. These include: 


e Ram pressure stripping: When a body moves with velocity vu, relative to a 
background gaseous medium of density p, the body experiences a ram pressure of 
magnitude ~ ov2,. This ram pressure strips gas from satellites and this gas mixes 
with the CGM of the central galaxy, contributing both mass and metals. In dense 
environments, ram pressure plays an important role in quenching star formation in 
satellite galaxies. This effect has been studied extensively in the context of galaxy 
clusters (e.g., and is theorized to produce “jellyfish” galaxies 
(e.g.,|Franchetto et al./2021). If the relative velocity of the satellite (or, better still, it 
full orbital history) is known, ram pressure can be exploited to infer the density of di- 
lute halo gas, as has been done using observations and modeling the Large Magellanic 
Cloud (LMC) around the Milky Way (Salem et al.[2015). 

e Tidal stripping: Tidal forces, which arise when gravitational forces are stronger on 
one side of a body than on the other, can also pull gas out of satellites. For a body 
of size AR, the tidal acceleration with which opposite parts of the body are pulled 
apart by the gravity of a point mass M at distance R is ~ GMAR/R?. We note 
that, because tidal forces scale as 1/ R?, tides between low-mass satellite galaxies that 
are near one another can be more important than tidal forces between a satellite and 
the central galaxy. For example, modeled the Magellanic Stream 
(a band of HI gas trailing the Magellanic Clouds) as due to LMC tides stripping gas 
from the Small Magellanic Cloud (other models suggest an important role for ram 
pressure in addition to tides, e.g. Lucchini et al.[2021). If 
the Milky Way were observed externally, the Magellanic Stream would appear as an 
important component of its CGM, so we must presume that some observed features 
of the CGM of other galaxies arise from similar tidal interactions. 
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e Satellite winds: Feedback in satellite galaxies can eject gas from the ISM of a 
satellite and into the CGM of the central galaxy. This effect is illustrated in the 
simulation shown in Figure where it is labeled “intergalactic transfer” because 
some of the gas ejected by satellites can later accrete onto the central galaxy. This 
transfer process can contribute up to ~ 1/3 of the baryons that end up as stars in 
Milky Way-mass galaxies in FIRE simulations (Anglés-Alcázar et al.[2017), although 
the importance of this mode of galaxy fueling differs in other simulations (bottom 
right panel in Fig. |7). Not all the gas ejected in satellite winds necessarily re- 
accretes onto galaxies. Winds from satellites can also affect the CGM by puffing up 
accreting filaments in which they are often embedded (see and by creating 
overdensities that promote the condensation of the cool gas via thermal instability in 
the CGM (Esmerian et al. [2021]. At a given time, the fraction of the total CGM mass 
contributed by winds from satellite galaxies can be substantial (e.g., up to ~ 20% 
inside the virial radius of L, galaxies in the FIRE simulations analyzed by Hafen et 
al. 2019). 


The different mechanisms listed above that remove gas from satellites are not mutually 
exclusive. The galaxy group containing M81 and M82 is a well known example of a system 
in which intra-halo, filamentary HI clouds are associated with strong galaxy interactions, 
and thus most likely involve tidal stripping (Chynoweth et al./2008). In this system, M82 is 
also well known for its prominent galactic A ee eei E kn Dog) so this 
an example where both gas ejection in a wind and tidal interactions shape the observed 
CGM. 

Even if gas mass losses by satellites are negligible, satellites can deposit into the diffuse 
CGM the gravitational potential energy they lose as they fall into, or orbit within, the halo. 
Ram pressure acts as an effective friction force and removes energy from the orbit, which can 
in principle go into heating the CGM (e.g., as wakes dissipate). Similarly, dynamical friction 
induces wakes that can dissipate in the CGM (El-Zant et al.|2004). 'These processes may 
contribute to the excitation of disturbances or turbulence in the CGM, and operate on gas 
clumps that accrete in the halo even if these clumps do not correspond to satellite galaxies. 
analyzed clumpy gas accretion in massive halos and argued based 
on simple estimates that “gravitational heating" by cosmological accretion delivered to the 
hot gas in the inner halo could potentially maintain star formation quenching in the long 
term. More recent high-resolution simulations that self-consistently include gravitational 
heating indicate that this process is not sufficient to quench galaxies (e.g., [Su et al.[2019). 
However, more work is needed to understand whether and where these processes may have 
interesting effects on the CGM. 


3. SMALL-SCALE PROCESSES: MULTI-PHASE GAS 


An outstanding problem in current large-scale simulations is that the amount of CGM cold 
gas is unconverged. It increases monotonically with resolution (Faucher-Giguére et al.|2016 


van de Voort et al.]2019] [Hummels et al.]2019] [Peeples ct al. [2019 


indicating that key physical processes remain unresolved (see Fig. |8). In this portion of the 


review, we will survey small scale processes in CGM gas. This frequently includes physics 
which is unresolved in galaxy or cosmological scale simulations, and is often the realm of 
idealized simulations or analytic theory. The list of relevant physical processes is vast, and 
similar to that in the ISM: magnetohydrodynamics, fluid instabilities, shocks, radiative 
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cooling, anisotropic conduction and viscosity, turbulence, cosmic rays, and gravity, to name 
a few. The more dilute nature of CGM plasma means that collisionality can be weak, and 
kinetic scale plasma processes can play a role. Entire textbooks could be devoted to some 
of these topics; we obviously cannot do them justice in a brief review. 

To focus our discussion, we lean on the striking abundance of atomic (T ~ 10* K) and 
sometimes even molecular (T ~ 10 — 100 K) gas in the CGM (Tumlinson et al.12017), even 
though virialized gas should be much hotter. Indeed, cold gas forms our main observational 
probe of the CGM: it is observed at high spectral resolution by quasar absorption line 
spectroscopy, and also in spatially resolved emission maps by Integral Field Spectrographs 
on large ground-based telecopes such as Keck and the VLT. Direct observations of the hot 
gas component, in X-ray and Sunyaev-Zeldovich measurements, are few and far between, 
except in hotter systems such as massive ellipticals, groups and clusters; dispersion measure 
(DM) measurements of fast radio bursts (FRBs; [Prochaska et al.[2019] ?) could eventually 
improve the situation. Given its observational prominence, and the likely importance of cold 
mode accretion in fueling star formation, we focus on physical processes relevant to cold gas 
in the CGM, and how it interacts with the hot phase. Parallels with the terrestrial water 
cycle are reflected in terminology (precipitation, condensation, evaporation). We consider: 


e Cold Gas Formation EES {3.2}. What is the origin of cold gas in the CGM? We 
have already discussed cooling flows which go through a sonic point (82.1.1) and 
cosmological accretion via ‘cold streams’ (2.1.4), both of which occur in limited 
halo and accretion rate regimes. Here, we will discuss three further possibilities: 
thermal instability of hot halo gas (‘precipitation’; (3.1.2), wholesale cooling of wind 
gas (83.1.3), and mixed-induced cooling condensation onto cold gas ‘seeds’ CES see 
below). 

e Cold Gas Survival and Growth CES Overdense cold gas cannot be supported hy- 
drostatically. It must either fall under gravity, or be flung out at high velocity. The 
resulting shear with hot gas should destroy the cloud via Rayleigh-Taylor and Kelvin- 


Helmholtz instabilities (Klein et al.|1994| Zhang et al.||2017). Yet, cold gas is seen 


in abundance outflowing in galactic winds, and also infalling as high-velocity clouds 
(HVCs; [Putman et al. [2012). We describe recent progress in understanding cold gas 
survival, and how it can even grow in mass. 

e Cold Gas Morphology (83.3). Small-scale structure (< 50 pc) in CGM cold gas has 
been inferred from photoionization modeling 


Stern et al.|2016 2019); it has been suggested that cold gas is a ‘mist’ 
(McCourt et al./2018). What is the topology of cold gas? Does it have a characteristic 


scale? Or is there structure on all scales? Like the IMF of self-gravitating clouds in 
the ISM, the mass function of pressure-confined cold gas in the CGM has fundamental 
observational and theoretical consequences. 

e Cold Gas Interactions: Turbulent Mixing Layers ($3.4). Phase boundaries are not 
infinitely sharp, but thickened by diffusive transport processes such as thermal con- 
duction, viscosity, and turbulence. The physics of these mixing layers are of great 
consequence: they govern the transport of mass, momentum and energy between 
phases. What sets these transport rates? These mixing layers are usually completely 
unresolved, even in idealized simulations. Does this preclude numerically converged 
transport rates? 

e Cold Gas Interactions: Cosmic Rays (83.5). Cosmic rays have energy densities com- 
parable to thermal gas in the ISM; an important role is entirely plausible in the CGM. 
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Figure 8: Cosmological simulations are unconverged in the CGM. Shown are column density 
projections of a simulated L, galaxy with different levels of spatial resolution. Enhanced 
halo resolution leads to enhanced neutral hydrogen content, with no sign of convergence. 


Adapted from (2019). OAAS. Reproduced with permission. 


'They can provide non-thermal pressure support, accelerate and heat gas; in recent 
years, their potential role in driving galactic winds has drawn significant attention. 
We briefly describe CR hydrodynamics, then draw attention to how small-scale cold 
gas can dramatically alter CR transport, and the spatial footprint of CR momentum 
and energy deposition. 


3.1. Cold Gas Formation 


3.1.1. Making Multi-Phase Gas: Local Thermal Instability. How does a multi-phase 
medium develop? The classic mechanism is thermal instability (Fiela]1965): slightly over- 
dense gas cools faster than its surroundings, loses pressure, undergoes compression and 
runaway cooling, until it reaches a new equilibrium. In the words of (1995), “This 
seems such an economical and elegant method to make cloudy media, one feels that nature 
would be inexcusably remiss not to have taken advantage of it at some point.” In stratified 
galactic halos, this is often dubbed ‘precipitation’, since the cool gas which forms rains down 
on the galaxy. In this review, we distinguish precipitation, discussed here, from ‘condensa- 
tion’, whereby hot gas mixes with pre-existing cold gas to form intermediate temperature 
gas, which subsequently cools (83.2). Both precipitation and condensation turn hot gas into 
cold gas via radiative cooling. T'hey can take place simultaneously: for instance, cold gas 
which forms via thermal instability (precipitation) will fall and shear against background 
hot gas, triggering condensation. However, their physics is quite different. Precipitation 
has recently been comprehensively reviewed by [Donahue & Voit] (2022), to which we refer 
the reader for more detailed discussion. 

Two points are worth noting. Firstly, local thermal instability typically presumes ther- 
mal equilibrium in the background medium; some form of heating is required. Otherwise, 
the entropy contrast between cool/hot gas often does not develop quickly enough, and a 
single phase cooling flow develops instead (2-1-1). Secondly, there is often some damping 
process which counteracts the fragmentation into a multi-phase medium. The dimensionless 
ratio tcoo!/tdamp, Where tcoo1 is the cooling time and taamp is the damping time, determines 
whether the medium is single phase (high tcooi/taamp) or multi-phase (low tcooi/tdamp)- 
Care must be taken to specify if Gool is the cooling time of hot, cool, or intermediate tem- 
perature gas; different timescales are relevant in different contexts. Examples of tgamp are 
the buoyancy time tpuoy in stratified environments (where tpuoy ~ £g; see below), and the 
eddy turnover time tturb in turbulent environments. 
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In a uniform medium in thermal equilibrium (£ = 0), the classic criterior{"| for thermal 


instability is (Field]1965): 


where £ is the net loss function (cooling minus heating per unit mass), S is entropy, and 
A is the variable held constant. The limits of when the sound crossing time is short (long) 
compared to the cooling time gives rise to isobaric (isochoric) thermal instability, where 
pressure (density) is held constant. For a cooling function A(T), with a local power law 
slope a = d(logA)/d(logT), equation (6) gives a < 2 (a < 0) as the criterion for isobaric 
(isochoric) thermal instability. For the range of temperatures relevant for thermal instability 
in the CGM (10° « T < 10'K), the gas is both isobarically and isochorically unstable 
(a < 0), in which case the isobaric mode has a shorter growth time, by a factor —yga/(2—a) 
(Field]1965), where yg = 5/3. There are also small ‘pockets’ in narrow temperature ranges 


where only the isobaric mode is unstable (Pfrommer|2013 2021). 


3.1.2. Precipitation: Local Thermal Instability in a Stratified Medium. In a gravitationally 
stratified medium, an overdense cooling blob will oscillate, due to buoyant restoring forces as 
it falls under gravity. The rapidly changing background a cooling blob experiences modifies 
thermal instability (Defouw[1970). showed that the medium should 
be thermally unstable if and only if it is convectively unstable, dlnS/dlnr « 0, where 
S(r) = P(r)/[o(r)]?/? is the radial entropy profile. Since virialized halos are convectively 
stable (at least in groups and clusters, where X-ray observations show entropy profiles which 
increase outward), this appears to rule out thermal instability in stratified halos. 

In fact, assumed that heating is a function of thermodynamic 
variables such as density and temperature (e.g., photoionization heating), but not an explicit 
function of position. Observations of AGN feedback in cluster cores means that spatial 
dependence is very likely; the same could be true of stellar feedback in galaxies. In that case, 


the|Balbus & Soker|(1989) criterion does not apply. In an influential set of papers, [McCourt] 
(2012) and (2012) showed that local thermal instability (commonly 


dubbed precipitation, since the cold gas then falls like raindrops) occurs when fcoo/tg < 
1,10 in planar (spherical) simulations; otherwise, the medium remains single-phase — see 
F igp] They also showed non-linear saturation at density amplitudes óp/p  (tcooi/te)~+ 

Cooling blobs are akin to driven damped oscillators, with the ratio of driving time (the hot 
gas cooling time tcoo1) and damping time (the buoyancy time tbuoy ~ (dlnS/ dlnr)- / ?ta) 
determining the existence and saturation amplitude of thermal instability. Damping arises 
due to non-linear g-mode couplings, which sap energy from the oscillating blob. There is 
also considerable observational evidence in groups and cluster] that a threshold value of 
min(teooi/t#) S; 10 is required for multi-phase gas (e.g., see summary in [Donahue & Voit] 
2022), and it has been hypothesized that massive galaxies and clusters naturally evolve, 
via feedback, to a state of marginal stability with respect this criteria (Voit et_al.]|2017). 


TA more SES criterion, [2(C/T)/8S|A > 0, drops the requirement of thermal equilibrium 
(Balbus 1986 

“For a enting view, see [McNamara et all McNamara et al. al.| ( 2016 ) 2016) and [Hogan et al. (2017), who argue from 
observational data that tcoo] alone, rather than t.oo1/tg@, determines the threshold for multi-phase 
gas; they contend that non-linear perturbations, such as uplift of cold gas from the central galaxy, 
seed cold gas condensation. 
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Figure 9: Nonlinear saturation of thermal instability depends only on the ratio tcooi/te 
in hydrodynamics. Each column shows images of the density in simulations with different 
values of the ratio Leed (fe: all panels show the density after several cooling times, when the 
simulations have reached a steady state. Thermal instability produces clouds of cold gas 


(with àp/p = 1) only when tooo € te. Adapted from (2012). 


These considerations are also important in the vicinity of accreting black holes, which can 
be potentially fed by precipitating cold gas (‘chaotic cold accretion’), rather than Bondi 
accretion of hot gas, which is inefficient and also predicts a cusp in hot gas temperature 
near the black hole, which is not seen [2018]. Cloud-cloud and cloud- 
torus collisions allow angular momentum loss. 

We have already seen the ratio teooi/t# before, which determines when halos complete 
virialization ({2.1). In the absence of heating, the competition between cooling and shock 
virialization of cosmological accretion determines global thermal stability, i.e. whether cool- 
ing flows (or stable virial shocks) develop. In the presence of heating, and background hy- 
drostatic/thermal equilibrium, the competition between cooling and buoyancy determines 
local thermal stability, i.e. fragmentation into a multi-phase medium. While both virializa- 
tion and buoyant oscillations share a common timescale tg, their physics is quite different. 

For precipitation, the ratio tcoo1/t# is more correctly written as tcool/tbuoy, Where thuoy ~ 
(dinS/dinr)~'/?tg. Typically (dinS/dinr) ~ O(1) in halo gas, which is strongly stratified. 
However, in environments with strong turbulent mixing, entropy cores can develop, and 
then there are no buoyant restoring forces: tbuoy — oo, and teool /t is irrelevant 
2017). While observational constraints are weak, flat hot gas entropy profiles often develop 
in simulated disk galaxies (Esmerian et al.[2021). Furthermore, global thermal equilibrium 
requires that any feedback loop which keeps the hot atmosphere stable must operate on 
timescales short compared to thermal timescales. This likely holds in the hot atmospheres 
of galaxy groups and clusters, where the AGN duty cycle (t ~ 10’yr) is shorter than, or 
comparable to the central cooling time[] However, the situation is likely much more dynamic 
and out of equilibrium in galactic CGM environments (e.g., [Fielding et al.[2017b). Thus, 
while the case for precipitation in massive halos (M > 1018 Mo) is strong, the extrapolation 
to lower masses is more uncertain, particularly given scant observational constraints on hot 
gas in such halos. In these environments, it is quite possible that buoyant restoring forces 
are secondary to other physics (e.g. turbulence). 

Indeed, the threshold value of Leed (fe can vary in different setuingd or have consid- 


9Many simulations adopt idealized heating where global thermal equilibrium is enforced by fiat: 
at each timestep, total cooling in each shell is calculated and an equivalent amount of uniform 


heating is added to that shell (McCourt et al.|2012| |Sharma et al.|2012). However, simulations 
which incorporate more realistic AGN feedback heating obtain similar results (Gaspari et al.12012 
Li & Bryan|2014a]. 


“It was initially thought that damping by buoyancy is geometry dependent, given different 


thresholds for planar and spherical simulations (McCourt et al.||2012 2012), but 
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erable dispersion (e.g., [Voit]2021). A particularly important variable is the amplitude of 
initial density perturbations, which is presumed small (0p/p « 1) in most studies. In a 
controlled set of simulations, showed that the threshold value of 
Lei (fo increases (i.e., thermal instability is more easily triggered) as óp/p increases; once 
óp/ p ~ O(1), then buoyant forces cannot quench thermal instability, and Leet (fe becomes 
irrelevant. Such non-linear density perturbations could arise from shocks, or the introduc- 
tion of low entropy gas into the halo via uplift from the halo center (e.g., cool gas dragged 
out by AGN blown bubbles), cosmological accretion, or winds from satellite galaxies. Such 
low entropy gas cools easily and seeds further cooling of hot gas, though in the context 


of this review, this falls under the rubric of ‘condensation’ (83.2) rather than ‘precipita- 
tion. Both and came to similar conclusions 
that tcooi/t# is a poor predictor of thermal stability when large amplitude perturbations 
are present, in TNG50 and FIRE simulations respectively. 

Other physics can also influence the development and saturation of thermal instabil- 
ity. For instance, magnetic fields damp buoyant oscillations (via magnetic tension) and 
destabilizes all scales below la ~ vatcoo! (where va = B/4/A4mp is the Alfven velocity) — 
i.e., thermal instability can potentially occur anywhere in the halo, independent of Leed (Te 
(Ji et al.|[2018). MHD effects appear for 6 = P;/Pg < 1000 (where P} is gas pressure 
and Pp is magnetic pressure), with B-fields enhancing the amplitude of thermal instability 
óp/ p c B -V/? even for very weak fields, surprisingly independent of field orientatior!!] By 
providing rotational support against gravity, and also deflecting the descent of a condensing 
blob via Coriolis forces, angular momentum can also suppress buoyant damping of ther- 
mal instability (Sobacchi & Sormanil/2019). 'Turbulence can both enhance and suppress 
thermal instability. Enhancement occurs via stronger density perturbations in compressive 
turbulencd??] as well as uplift of low entropy gas, which subsequently cools (Voit|2018). 
Suppression occurs via mixing of cold gas with hot gas; have argued 


that tcool/tturb S 1 is required for perturbations to condense, which we will see is indeed 
the criterion for multi-phase gas in a turbulent medium (see q3.4]. Cosmic rays are an 
interesting case. In simulations with horizontal fields perpendicular to gravity, they do not 
significantly change the threshold for thermal instability, although (by providing pressure 
support) they decrease the density, and increase the infall time of cold clouds 
2020). However, since field-aligned CRs cannot stream down their vertical gradient, there 
is no background CR transport or CR heating in the horizontal field case; these only arise 
as perturbations. For vertical fields with CR streaming, linearly unstable modes no longer 
oscillate like g-modes, but instead propagate at the Alfven speed 
[2020]. In the non-linear state, in some parameter regimes hydrostatic and thermal equi- 
librium can be maintained, but in others, a CR heated and driven wind develops — local 
thermal instability triggers a global wind instability (20222); see (3.5). 

Once dense cold gas is introduced into the halo, it cannot stay in equilibrium but either 
falls under gravity or is blown out by a wind. It therefore inevitably evolves as it shears 


subsequent work showed that this was an artifact of the how tcooi/te was defined in different 
setups, where teoo1/tg can either increase monotonically or display a minimum. The physics of the 
threshold is geometry independent (Choudhury & Sharma[2016| Meece et al.[2015). 
11 Horizontal field lines support overdense gas directly via magnetic tension, while vertical field 
lines confine over-pressured hot gas, which in turn supports overdense gas via pressure gradients. 
1?Even for subsonic turbulence, where p/p ~ M <1, where M is the Mach number, non-linear 
steepening can result in dp/p ~ O(1) and weak shocks. 
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against hot gas, potentially either growing (‘condensation’) or diminishing (destruction via 
hydrodynamic instabilities). Such processes are critical in regulating the cold gas content 
of halos. We discuss this in 


3.1.3. Global Thermal Instability in a Wind. Many galaxies do not have a static CGM, but 
exhibit outflows and inflows, where cold gas is observed. Without stable heating, global 
thermal instability can cool the entire CGM within a critical radius if the cooling flow 
goes through a sonic point ( , or beyond a critical radius in a wind (Wang) 1995} |Silich 


2004| [Thompson et al.|/2016 [Scannapieco|2017 Schneider et al.|/2018). 


et al.| (2016) noted that mass loaded winds always cool on large scales. This may seem 


counterintuitive, since in an initially adiabatic spherical wind (Chevalier & Clegg|1985), 


where v const, the density p x r ?, and cooling rates fall rapidly. However, as the 


wind cools adiabatically, it moves into a temperature regime (10? K < T < 10’K) where 


metal line cooling dominates, A c T~°". Since T x P/p x pt « r-*/5. we have 
teoo. X T/(nA) oe T [n x p 4/15. and since the advection time t x r/v œ r, this gives 
tesoi/ tay Eé p-19/15. which decreases rapidly with radius. [Thompson et al.|(2016) derived a 


‘cooling radius’ for the outflow, when Lei ~ taav, beyond which the flow becomes radiative: 


2.13 0.789 

walle a RLT? Qar 7 

Tcool ~ pe 2.92 10.3 ` ? d 
B? Msrr,10 


where a is the fraction of supernova energy thermalized in the hot plasma (Ewing = 
3x 104T erg s^ aMsrn), B= Mwina / MsFR is the mass loading factor, Ro.3 = R/(0.3kpc) is a 
characteristic radius for the injection region, Msrr.io = Msrr/(10 Moyr™!) is the star for- 
mation rate, and Qa, = 2/47 is the opening angle of the outflow. This model bypasses the 
need to entrain cold gas, which is born comoving, and produces cold gas in abundance, since 
essentially the entire wind cools to T ~ 10*K at red, From equation (7), it is clear that the 
viability of a cooling wind is sensitive to the kinetic injection factor o and mass loading fac- 
tor f: estimate that for B < Bmin © 0.6099 (Ro. al / MsrR,10)9?6, 
the wind remains adiabatic to large scales, since tcoo1/taav never falls below unity before 
T < 10?K, when the cooling function changes slope. Indeed, in high resolution simulations 
with high (8 ~ 0.5) mass loading rates, found wind cooling in good 
agreement with equation (7), but no wind cooling for simulations with 6 ~ 0.1 
[der et al.|/2020), consistent with above expectations. This sensitivity to mass loading also 
appears in semi-analytic models of multi-phase winds (Fielding & Bryan[2022). 


3.2. Cold Gas Survival and Growth 


3.2.1. Condensation Enables Entrainment and Growth. We now describe condensation, 
which like precipitation turns hot into cold gas. A combination of cold gas ‘seeds’ and 
velocity shear against the hot phase produces mixed intermediate temperature gas, which 
subsequently cools. The cooled gas retains its momentum, and in time the cold phase 
both grows in mass and comoves with the hot phase. This form of mixing-induced thermal 
instability thus enforces kinematic coupling between phases, and the high rate of conversion 
between phases has important implications for the baryon cycle. 

Condensation is closely related to the well-known 'cloud entrainment problem'. Atomic 
and molecular gas is observed outflowing (in galactic winds) or inflowing (as intermediate 
or high velocity clouds, IVCs/HVOSs), at velocities often comparable to virial velocities. 
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Figure 10: Clouds with toool,mix/tec < 1 can entrain, survive and grow in a wind. While the 
simulation on the left, where tcool,mix/tcc ~ 7.7 shows the expected behavior well known 
from previous studies, i.e., the destruction of the cloud on a few tec, the simulation to the 
right, where teool,mix/tcc ~ 0.077 shows that the cloud does not get destroyed but instead 
grows to a mass greatly exceeding its initial mass (m/s is the mass of gas with density at 
least 1/3 of the initial cloud density), and also gets entrained, with Av/vwina < 1. The 


grey disk in each panel illustrates the original size of the cloud. From|Gronke & Ob (2018). 


How can it withstand hydrodynamic instabilities? The problem is easily stated: if cold 
gas is overdense by a factor x = pc/pn, the timescale for acceleration by hydrodynamic 
ram pressure, taco ~ p/p ~ (pe R°?v)/ (pnv? R?) ~ xR/v, is longer than the timescale for the 


cloud to mix into its surroundings, the ‘cloud-crushing’ time, tec ~ tku ~ J/XR/v, by a 
factor tacc/tee ~ x1? ~ 10 — 30, independent of cloud size or velocity (Klein et al.| 1994 
Faucher-Giguere et al.|2012| Zhang et al.|2017). Here, pn, p. are the densities of hot and 


cold gas, R is the cloud radius, v the relative velocity between the cold and hot gas, and p is 
the momentum. Many simulations, both with and without radiative cooling, confirmed that 


clouds in a wind tunnel are destroyed before they are entrained (Klein et al.|1994 
2002| |Pittard et al.[2005 2009} |Scannapieco & Brüggen|2015 
& Robertson|2017). Here, we describe an emerging paradigm for cloud acceleration and 


survival where hot gas condenses onto cold gas, by mixing and cooling, thereby transferring 
mass and momentum to the cold phase. 

Cloud growth (rather than destruction) was first seen in wind tunnel simulations with 
radiative cooling by and [Armillotta et all (2016], in the context of 
HVC survival. (2018) found not only cloud growth, but cloud entrainment 
in a wind; they quantified a criterion for cloud survival, tooolmix S tec, which matched 
simulation results well. Here, tcool,mix is the cooling time of mixed gas at Tmix ~ V To Ty 


(Begelman & Fabian 1990 $3.4), where Te, Th are the cold and hot phase temperatures. 
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This criterion sets a lower bound on cloud size (Gronke & Oh/2018): 


5 


T Mwina x 
P3Amix,—21.4 100 


where Tu = (Ta/10* K), P = nT/(10? cm? K),  Amix,-21.4 = 
A(Tmix)/(10~?!4 ergem? s7!), Mwina is the Mach number of the wind, and we write 
Uwind = Cs,wind/VÍwind ^v Cs, cM winax T2, assuming isobaric conditions. In all cases where 
the cloud survives, it forms a cometary tail, similar to the observed head-tail morphology 
of HVCs (Putman et al. [2011], which grows in mass. Cloud growth and entrainment are 
intertwined: as hot gas condenses, it imparts its momentum, and accelerates the cloud. 
See Fig. [10] for an example. 


One can quantify growth and entrainment times (Gronke & Oh|2020a); analytic models 


match simulation results well. Interestingly, cloud growth does not stop once the cloud is 


Uwindtcool,mix 
R> : 


xi x 2pc 


entrained, even though there is little shear to drive mixing. Instead, growth is maximized for 
comoving clouds! In entrained clouds, mixing is driven by cooling-driven cloud pulsations, 
which arise for clouds which lose sonic contact as gas condensing onto the cloud cools, 
tcool «€ tse (where tse ~ Hie is the sound crossing time). The small loss of pressure balance 
causes the cloud to contract, overshoot and subsequently expand, which in turn drives more 
mixing and cooling, and subsequent pulsations. Cooling in a pressure confined cloud drives 
overstable acoustic oscillations governed by tcooi/tsc, just as cooling in a stratified medium 


drives overstable buoyant oscillations governed by ícoo/tm (83.1.2). Such pulsations are 
also seen in clouds in initially static surroundings which are perturbed ben die 
2020b). They are the linear, small-amplitude version of ‘shattering’ (83.3.2), where the 
entire cloud falls violently out of pressure balance due to cooling, creating highly non-linear, 
large amplitude oscillations which tear the cloud apart. Hot gas inflow velocities were found 


—1/4 


to be Umix ~ Cs,c(tcoot/tsc) , and thus growth happens on a timescale tgrow ~ m/m ~ 


Dc Ar / (p AUmix) ~ Xtsc. We dissect Umix, and the unintuitive scaling umix ox (teoo1/tsc) 4. 
in Surprisingly, mass growth is converged at fairly low resolution (~ 8 cells per 
cloud radius; [Kanjilal et al.[2021), even when the mixing surface area 
is manifestly unconverged, and thermal diffusion lengths or Field lengths are unresolved 
(indeed, even without explicit conduction); these issues are also discussed in 83.4 

The criterion teool,mix < tec, where Tmix ~ ew is an idealization which happens 
to match simulations with temperature floors of T ~ 10* K (due to photoionization) well. 
In reality, mixed gas has a broad range of temperatures. A more refined criterion based 
on the thermodynamic history of fluid elements, modeling how mixing and cooling modify 
entropy, would be more accurate, at the expense of greater complexity 
[2022). If gas cools to lower temperatures, different criteria (still all variants of a cooling 
length) potentially apply (Farber & Gronke| 7). The definition of cloud survival is 
also important. (2020) and (2020) argued that it is the cooling time 
of hot gas which matters, i.e. that tcool,hot < Occ (where a ~ 10, since it takes several 
cloud-crushing times to destroy a cloud) is required for survival. This is indeed true if one 
cares about the survival of original cloud material, as appropriate in some applications. 
However, the abundance of cold gas is not monotonic. After an initial decline, eventually 
the mixed gas cools, and the cold gas mass recovers and continues to grow[?] even if much 


of the original cloud material has not survived (Kanjilal et al.12021). 


13The cooling time of the hot gas cannot be the rate limiting step for cloud growth, since turning 
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The physics above has been demonstrated in 3D hydrodynamic simulations with tran- 
sonic, constant winds and x ~ 100 — 1000. Below, we discuss how other parameter choices 
and physics affect outcomes, but cover turbulence CES and cosmic rays (83.5) elsewhere. 
These are all issues of active research, and by no means settled. Of these, it is most ap- 
parent that infall under gravity changes survival criteria (83.2.5), and the impact of higher 
overdensities and wind Mach numbers is most uncertain (83.2.4). 


3.2.2. Condensation: Magnetic Fields. The hot medium B-field has a much larger impact 
than initial cloud B-fields (Li et al.]2020). Unfortunately, the plasma 8 of the CGM is still 
highly uncertain, although it is expected to lie between the ISM (8 ~ 1) and ICM (8 ~ 100). 
Effects are strongest when B-fields are GE to the direction of cloud motion, and 
‘magnetic draping’ (Lyutikov| {2006} [Dursi|/2007) takes place: the cloud sweeps up field 
lines, which are amplified to rough equipartition with ram pressure (€B,drape ~ OD wind Vind; 
where €B,drape is the magnetic energy density of the drape and a ~ 2). The amplified 
B-fields have two important dynamical effects. Firstly, magnetic drag couples the cold 


and hot gas, reducing acceleration times by a factor ~ (1 + 2/[f ina. M?]) ! compared 


to the hydrodynamic case (Dursi & Pfrommer|2008 2015). Secondly, B- 


fields suppress the Kelvin-Helmholtz (KH) instability via magnetic tension. This clearly 
suppresses mixing in both adiabatic and radiative (Ji et al.[2019) KH 
simulations. These two effects might be expected to have a major impact on cloud survival 
and mass growth rates. In fact, while cloud morphology changes dramatically, becoming 
much more filamentary, the survival criterion (equation EG and cloud mass growth rates 


appear relatively unchanged] (Gronke Gronke & Oh Oh Gronke & Oh|[2020a] [Li et all i et al. [2020] [Sparre et al.|[2020] 
[Jennings et al.|2022). For typical overdensities x ~ 100, magnetic drag only influences 
survival for strong background fields (wina S 1). Why does the strong draped field not 


suppress mixing completely? Although the B-field strength of the hot medium does not 
affect the energy density of the draping layer or the time it takes to grow (~ R/Uwina), it 
does affect its thickness, and thus the scale of the modes which are stabilized’?| For weak 
B-fields (high 8), only small scale modes are stabilized, and mixing can still occur, albeit at 
a reduced rate. Cloud entrainment still relies primarily on momentum transfer via cooling 
of mixed hot gas. 

Still, given that mixing is attenuated, the fact that mass growth rates are not strongly 
suppressed is puzzling; only for relatively strong fields (wina ~ 1) is significant suppression 
seen (Jennings et al.|[2022). One important clue is that growth rates are suppressed in 
MHD simulations of accreting cold streams, similar to KH simulations (N. Mandelker, B. 
Tan, private communication), and unlike clouds at identical Mach numbers. Evidently, 
clouds have additional degrees of freedom which change mass growth. This issue is still 


cooling on or off for T > 0.6 Th appears to have little impact on mass growth rates 
[Abruzzo et al.[2022]; emission is dominated by lower temperature gas. This suggests that 
mixing is responsible for the initial drop in hot gas entropy. 

4Grannow et al] claim significantly reduced mass growth rates in the MHD simulations, 
but they only followed cloud evolution for t Sj 3tcc. Their results are consistent with other MHD 


sims where growth rates rise at later times. 

15The draping layer has Le larape ^v R/(6aM3 ) and stabilizes modes A £$ S 10larape ~ 
R/M? ~ R/(BwinaM?2) , where Ma = v/va, Ms = v/cs are the Alfven and sonic 
Mach numbers. The total au Ster in the drape is E ~ €p drapeAldrape ~ pv? ARV, v? ~ 
€B,windV, ie. it is the magnetic energy of the wind displaced by the cloud. 
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unresolved, but two points are worth noting. Firstly, morphology is very different: clouds 
can grow extended tails and also fragment; the significantly larger surface area measured in 
simulations — due to lower cloud densities (see below) and filamentary morphology — could 
boost the mass accretion rate, even if mixing rates per unit area area reduced. Secondly, 
mixing in shear layers and cold streams is via the Kelvin-Helmholtz instability — if shear 
drops to zero, so does mass growth. By contrast, mass growth in clouds peaks when it 
is entrained, due to cooling-induced cloud pulsations. Thus, the suppression of the KH 
instability by B-fields may have a weaker effect. 

Note that plasma £8 and density can vary dramatically in growing magnetized clouds. 
Compressional B-field amplification of entrained wind material can lead to low density, 
magnetically supported gas [Nelson et al.|2020). Even though high 
density regions dominate by mass, these low density, low 6 regions dominate by area, and 
could account for unexpectedly low density cloud material inferred from COS observations 


(Werk et al.|2014). In a sense, non-thermal pressure support can accelerate simulation 


convergence, since it reduces the scales to which gas is compressed and needs to be resolved. 


3.2.3. Condensation: Thermal Conduction. Conduction does not significantly affect the 
structure of the cloud itself, since hot electrons can only penetrate a skin depth Askin ^ 
Aen / X ~ 0.037? «nc! pc, where Ae, is the electron mean free path in the hot medium 
(equation [14]. However, it can create a thick boundary layer of warm gas around the 
cloud, which affects mixing and cooling. One can show that in the cloud survival regime, 
T > Toit (equation [B] we should be in the classical diffusive (Spitzer) regime; for a hot- 
ter background medium, conduction becomes saturated (Li et at.|020). Also, the ratio 
Terit / ÀF ~ 0.5f- 1/2 MT; 9/4 (where Ar is the Field length, equation |15), so conduction 
need only be mildly suppressed for clouds which survive mixing to also evade evapora- 
tion; cloud crushing and evaporation times are comparable, for typical CGM parameters. 
Accordingly, simulations with isotropic Spitzer conduction show modified cloud-wind inter- 
actions [Armillotta et al.[2016), but not dramatic changes. 
Effects include cloud compression and slower entrainment (attributed to the smaller cross- 
section); whether conduction impedes or enhances destruction depends on Mach number 
and overdensity. 

Importantly, however, real conduction is anisotropic, since electrons gyrate around B- 
field lines; cross-field conduction is strongly suppressed. For this reason, simulations with 
field-aligned conduction find the effects of conduction to be much weaker; once magnetic 
fields drape over the cloud, conduction is strongly suppressed, and has little effect on cloud 


mass evolution, relative to MHD simulations (Li et _al.]|2020 2022). Con- 


duction can also be suppressed by modified electron scattering rates, e.g. due to whistler 
instabilities driven by heat flux 
[2021). Still, the influence of conduction on cloud growth is not fully mapped out, and there 
could still be surprises. For instance, in all simulations, conduction considerably increases 
the amount of warm gas in the simulation domain. It still arises in anisotropic conduc- 
tion, since downstream gas is not magnetically shielded by draping. Downstream gas does 
not influence cloud evolution in current setups with conduction. But if there are multiple 
clouds, or if the velocity field is turbulent and time-varying, the abundant warm gas will 
interact with cold gas, with stronger effects. 
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3.2.4. Condensation: Wind and Cloud Properties. Wind tunnel simulations are of course 

highly idealized. They frequently assume a spherical, warm (T ~ 10*K, x ~ 100) cloud in 

a constant, trans-sonic (M ~ 1) wind. What is the effect of relaxing these assumptions? 
Begin with cloud properties. A number of simulations have considered non-spherical 


clouds - e.g., with initial conditions extracted from a turbulent box (Schneider & Robertson 
2017||Liang & Remming|2020, |Banda-Barragan et al.|2019| |Gronke & Oh|2020a 


2020). These appear to have little effect, apart from somewhat faster growth or destruction 
in the cloud growth/destruction regimes respectively, presumably due to larger surface area. 


Dust-laden molecular gas is often observed outflowing at high velocity (Ménard et al.]2010 
2010 2014); the higher overdensitied!*] and different cooling 


regimes at lower temperatures are important modifications. Given their formation path- 
ways, dust and molecules are often thought originate from the host galaxy; if so, they must 
survive the mixing process without being destroyed by high temperature gaf] Simulations 
of the entrainment of molecular gas find a complex variety of outcomes, with cold gas sur- 
viving, destroyed, or transformed to T ~ 10^K gas (Farber & Gronke]2022). Importantly, 
entrainment was faster than the naive expectation tace ~ xr/v, due to a ‘cocoon’ of warm 
T ~ 10^K gas forming, which decreases effective overdensity. Indeed, molecular gas is al- 
most always surrounded by atomic gas; in galactic HVCs, the atomic gas mass can be larger 
by a factor ~ 10 (Lehner et al./2009). This would lead to acceleration more analogous to 
the ‘standard’ x ~ 100 case; the atomic gas potentially protects the dust and molecules 
from exposure to hotter gas. 

For highly supersonic winds, as expected in starburst galaxies, the cloud develops a 
strong bow shock; it therefore interacts with higher pressure post-shock gas. Simulations 
of high Mach number winds show tails and entrainment times which are longer by a factor 
~ (1+M) (Scannapieco & Brüggen|2015| 2022). This can be attributed 
to the compression by an oblique shock, so that x > x(1 + M), which also produces a 


streamwise pressure gradient (Scannapieco & Brüggen . One important outstanding 
issue is convergence: 3D simulations (Gronke & Oh|2020a) are not converged at high Mach 


Bustard & Gronke|2022). This is not 


understood, though it may be related to suppression of lateral instabilities in 2.5D. In the 


numbers; only 2.5D simulationd'*]are well-converged ( 


latter, growth times scale as tgrow cx M? for high Mach numbers, and the survival criterion 
Ícool,mix < tcc remains surprisingly robust (Bustard & Gronke[2022). 

'This points to a potentially important disconnect. The poster child for the galactic wind 
community is M82, where velocities are ~ 2000 km s! , the expanding wind evolves rapidly, 
and the Mach number increases with radius (since v const and the sound speed falls due 
to adiabatic cooling); for instance, M ~ 2 — 8 is possible, though the lower end of the 
range is likely relevant if most clouds are launched at small radii. Still, the very hot winds 
at the launch radius (T ~ 108K) also imply larger overdensities (x ~ 10*) than are often 
simulated. This contrasts with M ~ 1 — 2, x ~ 100 — 1000 in typical idealized simulations, 
which is more appropriate for relatively quiescent CGM conditiond"?] For strong starbursts, 


16 High overdensities x are challenging to simulate, as clouds typically develop tails of length ~ yr, 
so longer simulations boxes are needed. 
17Though molecule formation and dust regrowth is also potentially possible in wind-driven ra- 


diative shocks (Richings & Faucher-Giguére|2018a b). 
18Such simulations only have a limited number of cells and periodic boundary conditions in the 


‘0.5’ direction. 
19E.g., the infalling cold clouds simulated in Luminous Red Galaxy (LRG) hosts in TNG50, which 
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global simulations verify some aspects of the condensation picture, but do not agree in 
others. Using the GPU based Cholla hydrodynamics code 
for uniquely high resolution (~ 5 pc) global disk simulations of galaxy outflows over a large 
(20 kpc) domain, cold clouds were found to accelerate out to ~ 800kms ^! (Schneider et al. 
2020). They attributed this to mixing of hot gas, from the linear relationship between 
cool gas velocity and value of a passive scalar originally injected in wind material. As the 
value of this passive scalar in cold gas increases — indicating hot wind gas which has cooled 
out — the cold gas velocity rises. Specific cold gas momentum and energy fluxes increase 
monotonically with radius, as one might expect if cold gas is being entrained. At the same 
time, the cold gas mass flux is not monotonic — it rises then falls, contradicting expectations 
of steady cloud growth in adiabatically expanding winds (Gronke & Ohj/2020a), albeit for 
much lower overdensities and Mach number£??] An important clue is eae 
find cold gas to be under-pressured by a factor ~ 10 relative to the hot phase, unlike 
the pressure balance found in idealized simulations. While its origin is not understood-it 
could be an artifact of under-resolving cold clouds — underpressured clouds have different 
cooling times and are vulnerable to cloud-crushing shocks. More work needs to be done in 
the computationally challenging high overdensity and Mach number regime appropriate for 
starburst winds, particularly in understanding convergence. 


3.2.5. Condensation: Gravity. Besides outflowing in a wind, cold clouds also fall onto the 
central galaxy under gravity. Such clouds could originate from thermal instability in the 
halo, cosmological accretion, or galactic fountain recycling (Fraternali[2017). 'The question 
of whether such clouds can survive and grow is crucial for explaining observations of inter- 
mediate and high velocity clouds, as well as fueling of star formation in the disk 
et al.|2012). There is a crucial distinction between outflowing and infalling gas clouds. The 
latter gradually entrain, so destruction processes become weaker until the cloud comoves 
with the hot gas, and hydrodynamic instabilities are quenched. By contrast, infalling clouds 
never become comoving — they accelerate under the action of gravity, and have to survive 
unrelenting shear at the terminal velocity. HVC simulations often simulate a wind tunnel, 
which ignores this distinction, though some simulation include gravity (Heitsch & Putman| 
[2009], and find cloud growth [Grennow et al.|[2022). Condensation 
in the cloud's wake means that it is largely composed of cooled halo gas; the contami- 
nation increases linearly with time and from head to tail, impacting cloud metallicities 
(Heitsch et al.|2022). However, the survival criterion is more stringent (Tan et al.|2022): 
terow ~ XT / Umix < tec (rather than teool,mix < tec). Survival is only weakly sensitive to cloud 
size. Instead, it is sensitive to cooling; a critical pressure is required for growth and accretion 
induced braking. In our Galaxy, this corresponds to a distance ~ 10 kpc from the disk; above 
that, clouds should fall ballistically, and must be large enough to reach the 'zone of sur- 
vival’. The drag due to accretion-induced mixing Farag ~ mv balances gravity Fay ~ mg 
when UT grow ~ gm/Th ~ gtgrow, a terminal velocity which can be substantially smaller 
than traditional terminal velocities due to hydrodynamic ram pressure VT hydro ~ /2Rxg. 


Compared to the virial velocity wt ~ gt, we have UT, grow ^v Uvir(tgrow/t) < 1, i.e., infall 
is sub-virial (Tan et al.|2022). Sub-virial infall velocities are commonly observed in LRGs 


(Huang et al.|2016| |Zahedy et al.|2019) and galaxy clusters (Russell et al.||2016). These 


condense hot gas, grow, and survive over cosmological timescales (Nelson et al 20201. 


?0 Wind tunnels can mimic time-dependent winds by appropriately modified boundary conditions. 
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velocities are much lower than predicted terminal velocities from traditional ram-pressure 
drag, unless the ‘drop height’ of colds is fine-tuned (Lim et al.|/2008). They are in much 
better agreement with velocities from sotto nee (are To 

3.2.6. Condensation: Cold Streams. An important scenario for cold gas survival is su- 
personic cold streams in Maio > 101? Mg halos at z Z 2, which are potentially cru- 
cial for gas accretion onto central galaxies (82.1.4). They typically have Mach numbers 
M ~ 0.7 — 2.3, density contrasts x ~ 30 — 300 with respect to the hot gas, and stream 


radii ~ 0.03 — 0.3 times the virial radius (Mandelker et al.|2020a). Since they are poorly 


resolved in cosmological simulations, idealized simulations play an important role. Their 
geometry gives them properties intermediate between planar shearing layers (§3.4) and fi- 
nite sized clouds, supporting both surface and body modes. Surface modes are present 
only for M < Marit = (14- x 1/3)972? e.g., Merit © 1.3 for x = 100). 
Heuristically, pressure perturbations which drive the KH instability operate on a sonic time, 
which exceeds the flow time for supersonic flows. However, for slabs or cylinders of finite 
thickness, body modes — waves reverberating between stream boundaries — can destabilize 
the system, for M > 1. Thus, the system is always unstable, and adiabatic simulations 
find that streams disrupt and expand; thinner streams will not make it to the halo center 
Ee ith clon nag 
netic fields slow down but do not halt disruption (Berlok & Pfrommer 2019). However, in 
the presence of radiative cooling, streams which satisfy tcool,mix < ud equivalently 
which are larger than a critical radius r > rerit, survive and grow in mass, analogous to the 
cloud case (Mandelker ci al 2020a). 

Current idealized simulations assume an infinite stream in a uniform background; stream 
evolution as it encounters radially varying CGM gas has only been modeled analytically 
(Mandelker et al.|20205). More importantly, current setups include initial velocity shear 
but no gravity. Thus, streams decelerate as they grow. In reality, since they are falling 
under gravity, they should face unremitting shear at some terminal velocity, similar to 
infalling clouds ({3.2.5). Indeed, unlike idealized simulations, cosmological simulations find 
that streams reach a roughly constant terminal velocity, in strong disagreement with free- 
fall velocity profiles [Goerdt & Ceverino|2015), although the role of 
numerical viscosity needs to be carefully assessed. The lack of deceleration could change 
the survival criterion for streams, as it did for clouds. 


3.3. Cold Gas Morphology 


3.3.1. Why Does Cold Gas Morphology Matter?. Is there a characteristic scale for cold gas, 
and why should we care? As we have already seen, size affects the survival and entrainment 
of cold gas; only clouds with r > ferit © Utcool,mix/ JX are able to withstand hydrodynamic 
instabilities in a wind. Cloud size also dictates mass growth rates, and the impact of thermal 
conduction CES It affects the interpretation of a host of observables, ranging from column 
densities, kinematics, cold gas mass fractions, and radiative transfer, in large part because 
cloud size sets the ratio of volume to surface area. For instance, cloud size/column densities 


?lHere, tshear is defined to be the time in which the mixing layer expands to the width of the 
stream in the non-radiative case. For CGM parameters, this typically gives timescales comparable 
to the cloud-crushing time tec ~ x1? Rfv, though in general it differs. 
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Cold Clouds Form Small Initially. . . Even big clouds 
break apart 
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Figure 11: Shattering via thermal instability and mixing. Left: An initial perturbation 
cooling down to 10*K. When the domain size is large (~kpc), the perturbation shatters 
into much smaller fragments with a characteristic scale cstcoo1 ~ 0.1 pc. Middle: a ~ 10* K 
cold cloud moving through an ambient medium at ~ 10’ K. The bottom shows a zoom-in 
of the red box in the top panel; even this apparently empty region of the wake is full of 


tiny cloudlets of cold gas. Adapted from (2018). Right: Shattering on large 


scales. The top shows temperature map of merging halos at z ~ 4. The strong merger shock 
leads to thermal instabilities in the post-shock region, as in the rectangular box. Zooming 
in on this box (bottom) shows that cold gas (T ~ 2 x 10* K) in the sheet has shattered 
into ~kpc scale fragments, embedded in T ~ 10° K gas. Adapted from 


(2019b). (C'AAS. Reproduced with permission. 


affect ionizing and Lya photon escape (see {2.1.4] for a discussion of spatially extended Lya 
emission). The interpretation of kinematic line profiles changes if absorption comes from 
discrete ‘clouds’, or a ‘fog’ of droplets entrained in hot gas: non-thermal broadening would 
either trace small-scale cold gas turbulence or large-scale hot gas kinematics respectively. 
On a pragmatic note, a characteristic scale potentially sets the resolution at which numerical 
simulations should converge. 

In recent years, growing observational evidence for small-scale structure in CGM cold 
gas has emerged. From photoionization modeling, maximum cloud sizes | ~ Ny/n ~ 35 pc 


are inferred from observations of the CGM of z ~ 2 — 3 galaxies 2015 
2016); the cloud size could be smaller if the observed column is the result of 


intersecting many (fa > 1) cloudlets along the line of sight. Another argument comes from 
the surprisingly large area covering factor fa ~ O(1) of CGM cold gas inferred from both 
absorption and emission line studies, despite the high overdensities (x ~ 10?— 10°; [Hennawi] 
[et a1.[2015). Even if the mass fraction of cold gas is large fm ~ O(1), its volume fraction 
fv ~ fu/x ~ 107° is tiny, and it is hard to explain its apparent abundance, unless it is 


somehow arranged in a thin shell. (2018) argued that just like terrestrial 
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fog, this becomes explicable if cold gas is widely dispersed in tiny droplet4??] if there are N 
droplets of size | in a halo of radius R, fA ~ NI’ /R?, fy ~ NDP/R3, so fA/ fv ~ Bil zl 
A fog can also explain large, highly suprathermal line widths (e.g., o ~ 1000kms~' in CII 
absorption; [Hennawi et al./2015). If this reflects turbulence within the cloud, they should 
be torn apart, but for a ‘mist’, it just reflects background hot gas motions. Very small scale 
structure is also observed in HVCs, quasar BLR, BAL regions (see references in [McCourt] 
let_al.|[2018), and also in the ISM (Stanimirović & Zweibel|2018]. A ‘fog’ could explain 
the surprising success of the uniform slab model in explaining Lya spectra from galaxies 
Radio scintillation could potentially probe the ‘mist’ in the 


CGM (Vedantham & Phinney 


the ISM (Armstrong et al. 


, just as it probes small-scale density fluctuations in 


3.3.2. Fragmentation and Coagulation: Physical Processes. How can gas fragment to small 
scales? There are at least two possible mechanisms: (i) if the cooling time Le) falls far 
below the sound-crossing time tsc ~ R/cs in a cooling cloud, causing it to become strongly 
underpressured relative to surrounding hot gas, Pueua « Phot, the cloud-crushing shock 
can ‘shatter’ the cloud (McCourt et al.|2018]. 'These authors argued that similar to the 
Jeans instability, where gravitational fragmentation leads to a characteristic scale Aj ~ 
cst, fragmentation in cooling, pressure-confined clouds imprints a scale Leloudlet ~ Cstcool ^v 
0.1 pc (n/cm ?)^!, where the quantity cstcoo is evaluated at its minimum at T ~ 10*K. 
While this length scale depends on the ambient pressure confining the clouds, the column 


density through an individual fragment Neloudlet ~ 10!" cm~? 


is essentially independent of 
environment. Cooling of large (R >> cstcoo1) clouds which subsequently ‘shatter’ is easily 
triggered by large non-linear perturbations (e.g. via turbulence). Shocks are particularly 
effective at ‘shattering’ cold gas, either via compression of pre-existing cold gas 


[et a1.[2002), or fragmentation in the radiative shock (Mandelker et al.|2019b). Shattering 


can also arise during initially isobaric linear thermal instabilit once the cloud loses sonic 


contact with its surroundings as it cools (Gronke E On[2020b) [Das et al.[2021). (ii) Another 
avenue for producing small cloudlets is via mixing in Kelvin-Helmholtz instabilities; the tails 
of clouds in winds show a plethora of dense small clumps 
2019), which have a column density distribution which peaks at 
Notoudtet ~ 10! em? (Liang & Remming]|2020), corresponding to the characteristic scale 


l ^ cstcooy suggested for cooling driven fragmentation. Figure [11] shows fragmentation by 
these two respective processes. ‘Shattering’ due to thermal instability in halo merger shocks 


can produce very low metallicity Lyman-limit systems, far from any galaxy (Mandelker et al. 


2019b). Numerical resolution obviously affects the onset and efficiency of shattering, but it 


can even affect the multi-phase nature of gas. |Mandelker et al.| (2021) find that if cstcooi at 


22These arguments only apply to high overdensity gas, as suggested by photoionization modeling. 
For low overdensity gas, as reported in z ~ 0.5 COS observations (Werk et al.|2014), the area 
covering factor of large, low density clouds can be substantial. While these low densities — about 
an order of magnitude lower than expected — is still a puzzle, non-thermal pressure support could 
play a role. Indeed, such low densities and high covering fractions are seen in simulations where 
magnetic support dominates (Nelson et al.12020). 

?3By contrast, and somewhat counterintuitively, linearly unstable isochoric modes do not appear 
to break up into small pieces. In 1D simulations, such clouds are instead compressed to high density, 


and oscillate until they regain pressure balance (Waters & Proga|2019b 2021), though 


this needs to be verified in 3D simulations. 
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T ~ 10°K is unresolved, gas ‘piles up’ at T ~ 10°K as further cooling becomes inefficient. 

In detail, there are subtleties. In idealized simulations, ‘shattering’ does not appear to 
take the form of hierarchical fragmentation, as in the Jeans instability. Instead, the cloud 
is strongly compressed by its surroundings, overshoots, then ‘explodes’ into small pieces 
in a rarefaction wave [Gronke & ORO). 
For break up, the density perturbation must be highly non-linear (6p/p >> 1) before losing 
sonic contact with its surroundings. Strong density inhomogeneities, attributed to Raleigh- 
Taylor instabilities or small scale isobaric thermal instability?4] 
et al.|2021) are crucial to fragmentation: an overpressured, uniform cloud does not fragment. 
Crucially, break-up is only seen when the final overdensity of the cloud is xr Z; 300, otherwise 
the cloud recoagulates (Gronke E On|2020P). While a density threshold is intuitive, the 
value xr Z; 300 is not understood from first principles. If this high density threshold is 
robust, break-up via cloud crushing instabilities only happens in very hot winds and ICM- 
like conditions (Ti ~ 107K, x ~ 1000), rather than the CGM (Th ~ 10°K, x ~ 100), 
unless background turbulence can disperse droplets before they recoagulate. As for the 
production of cloudlets via mixing, this is only important for clouds in the destruction 
regime. There, given enough time, the entire cold gas mass (including the cloudlets) would 
mix into the hot gas. For larger clouds in the entrainment/survival regime, after an initial 
period of fragmentation, cloudlets in laminar winds are ‘focused’ onto the cometary tail, 
and the cloud remains monolithic. Indeed, there appears to be a fundamental contradiction 
between ‘clouds’ (which have to be large, r > ferit ~ 10cstcoot to survive) and a ‘mist’ of 
tiny clouds. In simulations, cold gas in a mist only survives by coalescing to form larger 
clouds. 

Thus, one cannot neglect coagulation, which competes with fragmentation. There are at 
least two kinds of coagulation. The most obvious is direct collisions, whereby two colliding 
cloudlets stick to one another, similar to dust grain growth in protoplanetary disks. It is 
heavily influenced by turbulence. While there is an enormous literature on turbulent coag- 
ulation (beginning with the seminal paper by [Saffman & Turner[1956], given its relevance 
to terrestrial cloud formation, there have not been studies for CGM-like conditions; more 
work is needed. A second mechanism is coagulation due to the advective flow generated by 


hot gas condensing onto a cold cloud (Elphick et al.11991/|11992| |Koyama & Inutsuka|2004 
Waters & Proga|2019a). The inflow velocities in these studies was set by thermal conduction 


(or numerical diffusion). If inflow is instead set by hot gas cooling in turbulent mixing lay- 
ers, from mass conservation the inflow velocity is vin ~ vmix(r/ra) ^, where a z 0, 1,2 for 
plane-parallel (e.g., a semi-infinite turbulent mixing layer), cylindrical (e.g., the cometary 
tail behind a cloud), and spherical geometry respectively. Cloudlets feel a gentle ‘breeze’ 
from this advective flow, and entrain, just as in standard wind tunnels (Gronke & On[2022). 
At first blush, it seems this should be utterly negligible: since Vmix ~ Ce, the corresponding 
Mach number is M ~ (umix/Cs,h) (ra/d)? S 1072, which should be easily overwhelmed by 
turbulence. Yet, it is undeniably observed in simulations: cloud debris coagulates onto a 
cometary tail in a wind, and small droplets coagulate into clouds. The geometric dilution 


2{Das et all argue from 1D simulations that clouds must be isobarically unstable but 
isochorically stable for fragmentation (this only happens in narrow temperature ranges; see [31.1]. 
Otherwise, small scale isobaric modes within the cloud are suppressed by compression of the un- 
derpressured cloud. This interesting hypothesis needs to be tested in 3D. They found this to be 
true only for linear density perturbations; the distinction between isochoric and isobaric stability is 
irrelevant for non-linear perturbations. 
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can be offset by a large increase in surface are?" and when the mass fraction of cold gas 
is high, increasing inertia, the local rms velocity dispersion can be comparable to c.c. 


3.3.3. Multi-Phase Turbulence: Structure on all Scales?. Individual clouds, and ‘mist’ in 
particular, do not have to be long-lived; they can be continuously created and destroyed, 
if there is a long-lasting supply of cold gas. This appears to be what happens in turbulent 
flows, where large growing clouds continually shed small clouds, which then mix into the 
hot medium (Gronke et al./2022). Crudely, cold gas condenses at large scales, which then 
cascades to small scales. The survival criterion r > ferit is similar to wind tunnel simula- 
tions, although turbulence makes survival close to rerit highly stochastic; only clouds with 
T > Terit are assured of survival. The velocity structure function of cold and hot gas are 
similar, which implies cold gas can be used to trace hot gas kinematics. The mass spec- 
trum is a scale-free power law, dn/dm œ m ?, with roughly equal mass per logarithmic 
interval; similar power-law tails are also seen in other contexts (Li & Bryan|/2014b). The 
entire system grows exponentially, on a timescale tgrow ~ \/(L/Vturb)tcoo! — similar to that 
of turbulent mixing layers (q3.4] here, L is the stirring length and vturb is the rms turbulent 
velocity), for good reason: — it is a TML writ large, with very similar physics, including 
self-similarity. While the small-scale ‘mist’ and large-scale ‘clouds’ have comparable mass, 
‘clouds’ dominate mass growth and ‘mist’ dominates areal covering fractions — and hence 
are much more likely to be observed along a random line of sight, in absorption line spec- 
troscopy. Viewed in this light, the question of whether cold gas is ‘cloudy’ or ‘foggy’ is 
misguided. It is both; there is structure on all scales, and the most relevant scale depends 
on properties we probe (column density, mass, emission measure, covering fraction, etc). If 
the mass distribution is a roughly scale-free power-law, as in the stellar IMF, cut-offs (such 
as Terit OF Cstcoo!) can still reflect important physical processes. Just as for the IMF, the 
cloud-size distribution is set by a complex interplay between accretion, fragmentation and 
coalescence, except clouds are pressure confined rather than gravitationally confined. The 
origin of the IMF is still debated in the star formation community, and is by no means a 
settled question. Given that example, and weak observational constraints on cloud sizes, we 
likely have a long way to go towards a robust theory of cloud masses in the CGM community. 
For example, magnetic fields could substantially change cloud morphology and the mass 
PDF, as they do in laminar flows (83.2.2). The interaction between turbulence, shocks and 


gravitational infall can produce a complex, rain-like morphology (Banda-Barragán et al. 
2021). 


3.4. Cold Gas Interactions: Turbulent Mixing Layers 


3.4.1. Physics of TMLs. Gas phases exchange mass, momentum and energy at phase bound- 
aries, via microscopic (thermal conduction, viscosity) and macroscopic (turbulent) diffusive 
processes. The mean free path of turbulent eddies is much larger than electron or ion mean 
free paths. Thus, turbulent diffusion tends to be more efficient, as demonstrated in its role 
in stellar heat transport (turbulent convection), or momentum transport (turbulent viscos- 
ity) in accretion disks. At phase boundaries where there is velocity shear or differential 


?5'lT here is an interesting analogy between the coagulation force, Fig ~ pv2 A1 A2/ (4nr?), and 
gravitational forces Fio ~ GA) M2/r? (Gronke & Oh [2022). While both have a 1/r? scaling, the 
former can be enhanced by breaking up into small droplets, which greatly increases the surface area 
and speeds up coagulation. By contrast, mass is conserved under fragmentation. 
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Figure 12: Some properties of turbulent mixing layers (TMLs). (a) TMLs have abundant 
high ions, such as OVI. Adapted from (2010). OAAS. Reproduced with 
permission. (b) The Damkóhler number (equation|9) determines if the TML is single-phase 

(Da «1; weak 2. C regime) or multi-phase (Da> 1; strong cooling regime). Adapted from 
Ge (2021). ( ) In the fast cooling regime, the TML phase boundary has a highly 
wrinkled, fractal oe Adapted from [Fielding et a) (2020) . OAAS. Reproduced with 


permission. 


acceleration, turbulence is driven by Kelvin-Helmholtz or Rayleigh-Taylor instabilities. In 
recent years, significant progress has been made in understanding turbulent mixing layers 
(TMLs), particularly the rich interplay between turbulence and radiative cooling, usually 


in idealized planar shearing setups (Kwak & Shelton|2010 2011 
2012 2015] [Ji et al.|[2018] [Fielding et al.[2020] [Tan et al.|2021 


[Chen et al.[2022). 'This physics lies at the heart of condensation, which enables cold 
gas survival and entrainment (83.2). It has close parallels with turbulent combustion fronts 
Se e al [T98). 

In the TML, hot and cold gas mix to form intermediate temperature gas, which sub- 
sequently cools. The net result is that hot gas is converted to cold gas. An early analytic 
paper suggested that TMLs are characterized by a temperaturq?*] Tmix ~ WEN ER ? and 
width | ~ vitcoot,mix, where vs is the turbulent velocity and tcool,mix is the cooling time of 
mixed gas (Begelman & Fabian]1990). Despite its crude formulation, tcoo1,mix appears to 
determine thresholds for cold cloud survival (equation[B) and governs the boundary between 
single and multi-phase TMLs (equation [). Suprisingly, however, the commonsense 'cooling 
length’ prescription | ~ vtc;oo;, which works well for radiative shocks, and predicts column 
densities N ~ nlcoo1 which do not depend on density (since fcoo1 « 1/n) turns out to be 
incorrect for TMLs (Ji et a1. [2019]. 

A key quantity is the inflow velocity umix of hot gas into the mixing layer. This deter- 
mines the mass flux m ~ Op Gent, and enthalpy flux é (5/2) Pu; [1 + M?] into the TML. 


?6 We can derive this result from a slightly different viewpoint. Suppose the hot and cold gas have 
a turbulent velocities vh, vc respectively. The cold and hot phases share kinetic energy, get? ~ pni, 
so that ve ~ vy / Xx", where the density contrast x = pc/pn. TMLs are isobaric, Py ~ Pe. This 
implies that whenever different phases mix, the hot phase dominates the enthalpy flux: Ej / Ec 
Py vy / Pcuc ~ xi? >> 1. Conversely, the cold phase dominates the mass flux rnp rie ~ py vy / pcvc ~ 
x V 2 « 1. These ratios have been confirmed in numerical simulations GE If hot 
gas provides enthalpy and cold gas provides heat capacity, mixing results in a mean temperature 
T ~ rhy de ~ pup d pete ~ (ThTe)*/?. 
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The factor M? represents kinetic energy of the hot gas, which is eventually thermalized. 
Since the inflowing hot gas mixes and cools, in steady state m gives the cold gas mass 
growth rate, and é gives the TML surface brightness. Recent work has found the following 
Umix properties (some of which are shown in Fig [12): 


e TMLs are roughly isobaric. Pressure fluctuations are small, and tend to decline with 
increasing resolution. 
e The inflow velocity is of order the cold gas sound speed Umix ~ cs,c. Thus, for Te ~ 104 


K, to order of magnitude vui. ~ 10km s 1. 


This may seem surprising, since this 
velocity is low compared to typical shear velocities, vn ~ 100 — 1000kms~!. The 
fundamental reason is that the inertia of the cold gas limits the rate at which the 
hot/cold interface deforms and the two phases mix. For a Kelvin-Helmholtz time 
tku ~ J/Xl/vn, the mixing velocity is mix ~ l/tn ~ vn/A/X ~ JMncs,c. Of course, in 
detail, Vmix must also depend on cooling rates, as we discuss below. 


TMt, multi-phase) or 


e The inflow velocity has different scalings in the fast (Umix ox t.i; 


slow (Umix X gue. single phase) cooling regimes. Just as tcoo1/t# determines whether 


a stratified medium is single or multi-phase, the Damkóhler number (Damkoóhler 


1940): 


t L 
Da- turb _ ; 9. 
tcool Ut,ctcool (Tmix) 


where L is the outer scale of turbulence, determines if the TML is single (Da< 1) 
or multi-phase (Da> 1) (Tan et al.]2021). Turbulent heat diffusion can be charac- 
terized by a diffusion coefficient Diu» ^ Del, For a single phase TML, umix ~ 
d T eet ex s (here, we abbreviate tcool,mix with tcoo1). On the other hand, 
when the TML is multi-phase, only a small fraction of the volume is actively cooling. 
Most gas is either in the hot or cold stable phase; only their relative fractions change 
with depth in the TML. The effective cooling time is the geometric mes of the 
cooling time and eddy turnover time Zoo ~ CORE NE) nd 2. the expected scalings 


in emissivity € ~ P/Tcool X (vc ite? ? have been directly verified in numerical 


simulations a ). In this strong cooling case, which is most relevant to 
GL 


condensation 
1/2 4 
Deurb / 3/4 L d 
Umix ^v - ^ UC. " 10. 
Tcool ` Leet 


Since they scale with Vmix, the mass growth rates and surface brightness follow the 


—1/2 ,—1/4 
same loool , Lat 


scalings. The fast cooling regime is the main regime of interest 
for us, since the criterion for cloud survival tec/tcool,mix > 1 (83.2) is equivalent to 
Damix > 1. 

e Convergence in net cooling and mass transfer rates only requires resolving large-scale 
motions. In the multi-phase strong cooling regime, the interface is highly wrin- 
kled, greatly increasing the turbulent front surface area Ar from its laminar value 


Ar. |Fielding et al.| (2020) showed from their simulations that it is a fractal, with 


27This is a common random walk result: the effective optical depth for a photon in a medium 
with scattering and absorption optical depths Ts, TA respectively is T« ~ J/TSTA. Another example 
is the Field length Ap ~ VAecvetcoo1, which is the effective mean free path of a thermal electron. 
It is set by the geometric mean of the electron's elastic (~ Ae, set by Coulomb interactions) and 


inelastic (~ fe Lee), where ve ~ J/3kT'/me is the thermal velocity) mean free paths. 
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Ar/ Ar, = (A/L)?-P, where A is the smoothing scale and D = 2.5 is the fractal di- 
mension measured in their simulations. Since the area (and volume) of the cooling 
region is resolution dependent, one might expect the total cooling rate to be resolu- 
tion dependent. However, simulation surface brightness and mass entrainment rates 
converge at surprisingly low resolution 
[2021], even when the cooling length cstcool is unresolved, and without explicit thermal 
conductio when the cooling interface is one cell thick. Note that for many prac- 
tical applications (e.g., cloud entrainment, SES convergence in these quantities is 
equivalent to convergence in mass, momentum and energy transfer between phases. 

As[Tan et al (2021) argue, we are already familiar with the resolution of this conum- 
drum. When we stir cream into coffee, the mixing rate is independent of molecular 
diffusion rates, even though the latter is ultimately responsible for fine-grained mixing. 
In Kolmogorov turbulence, turbulent dissipation rates e zz pv?/I are independent of 
viscosity. Microscopic diffusion (or viscous) times are extremely long at macroscopic 
scales. Instead, the cream (or kinetic energy) cascades from large to small scales, until 
the timescale for microscopic processes becomes short enough to take over. Similarly, 
turbulence increases the area of the fractal cold/hot interface until (numerical or 
physical) thermal diffusion takes over. However, the mixing rate, and thus overall 
cooling rate, is independent of the details of the fractal surface. Instead, as in the 
coffee cup, the rate of mixing is set by the eddy turnover time at the outer scale. This 
is the scale which needs to be resolved, which is relatively undemanding. Thus, for 
instance, cloud growth and entrainment is converged for clouds which are resolved by 


~ 8 cells; the TML is completely unresolved (Gronke & Oh|2018). 


These properties assume M S 1, hydrodynamic TMLs. Supersonic TMLs have been 
studied in the slow cooling regime (Yang & Ji[2022]. 'The TML separates into two zones: 
a Mach number independent zone (similar to what we have discussed), plus an expanding 
turbulent zone with large velocity dispersion. Turbulent dissipation dominates over enthalpy 
advection, and reverses the sign of mass flux: cold gas evaporates. In MHD simulations, 
B-fields amplify to quasi-equipartition with turbulence, at which point magnetic tension 
suppresses mixing (Ji et al.[2019). These effects potentially strongly curtail cold gas growth. 
However, due to the influence of finite size effects, they may be better examined in wind 
tunnel simulations of a macroscopic cloud or filament, rather than planar mixing layers. 
In supersonic flows, a standoff bow shock forms, and the cloud interacts with subsonic 
postshock gas (q3.2.4). Wind tunnel simulations with B-fields show cloud growth which is 
surprisingly similar between MHD and hydro simulations, unlike TML simulations (93.2.2). 


3.4.2. TMLs: Observational Predictions. A potential observational diagnostic of TMLs are 
high ions like OVI, which peak at T ~ 10°K in collisionally ionized gas. Such gas should have 
short cooling times (~ 10 Myr in CGM conditions), but is Ee seen in abundance 


in QSO sighttines (Tamron st aL UT Prochasia star) and HVCs 


[2014]. High ions naturally arise in steady state T'MLs. Moreover, the surprising un 


a velocity centroids in low and high ions 


2021) arises naturally in TMLs, since different portions of a mixing layer participate 


?8]f thermal conduction is implemented, it does not affect cooling rates and Vmix, as long as 


Dart ~ vt,cL > Deonauct = &/(pcp) (Tan et al.l 2021). 
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in the same large scale turbulence. 

However, predicted column densities of individual TMLs are significantly lower than 
observations, typically by two orders of magnitude (Slavin et al. 
Ji et al.|2019). Non-equilibrium ionization (Kwak & Shelton|2010| Ji et al./2019) and 
photoionization (Ji et al.[2019) can only increase column densities by moderate amounts. 
Thus, if TMLs are responsible for observations, sightlines must pierce hundreds of mixing 
layers— which may be plausible if the CGM exists as a ‘fog’ of tiny cloudlets (see $8.3). 

Another frequently used discriminant between theoretical models (conductive interfaces, 
cooling flows, etc; e.g., [Wakker et al.[2012) is observed lines ratios such as N(SiIV)/N(CIV) 
or N(CIV)/N(OVI). It is important to realize that line ratios are not uniquely predicted 
by TML simulations. Instead, they arise from the temperature dependence of thermal con- 
duction, &(T). Although macroscopic quantities like mass, momentum and energy transfer 
rates between phases are independent of thermal conduction, line ratios depend explicitly on 
the temperature PDF P(T) of the TML, which depends on small-scale interface structure 
and hence does depend on thermal conduction (Tan & Oh[2021). Neglecting thermal con- 
duction is equivalent to assuming temperature-independent thermal (numerical) diffusion. 
A simple 1D semi-analytic model for conduction fronts reproduces simulated temperature 


PDFs and line ratios for different conduction laws remarkably welf?] (Tan & Oh|2021). 


3.5. Cold Gas Interactions: Cosmic Rays 


Although they only represent a billionth of all particles, cosmic rays have an energy density 
comparable to turbulence, thermal gas and magnetic fields in the Milky Way ISM 
et al.|2015). They are accelerated primarily at supernova shocks, where ~ 10% of supernova 
kinetic energy is converted to cosmic rays, although secondary acceleration in the halo via 
termination shocks (Jokipii & Morfill|1987] Bustard et al.|2017) 
or turbulence (Drury & Strong|2017| Bustard & Oh|2022) is also possible. Since Coulomb 
interaction timescales are long, radiative losses by electrons do not affect CR protons. Thus, 
CRs are unaffected by radiative cooling. CR scatter orders of magnitudes more than photons 
(e.g., ~ 10° times before leaving our Galaxy), enforcing tight coupling and making them 
a better candidate for mediating feedback (Socrates et al.| ). Although the possibility 
of CR driven winds was noted early on (Ipavich|1975| |Breitschwerdt et al.|1991), in recent 
years there has been a veritable explosion of work on CR driven winds (e.g., 
[Quataert et al./2022b). If the halo becomes CR dominated, the pressure 


support provided by CRs can potentially allow the cool, photoionized (T ~ 10^K) phase to 
become volume-filling, with abundant OVI (Ji et al.|2020), and suppression of virial shocks 
(Ji et al.|2021). Even if the gas remains multi-phase, CR pressure support in cold clouds 
reduces their density and increases buoyancy, increasing free-fall times 
and altering kinematic absorption-line signatures (Butsky et al.[2022). Heating by CRs can 


also offset radiative cooling, though this has mostly been modeled in the ICM (Guo & Oh 


aco rommer uszkowski et al. ). Regrettably, this body of work is 
2008| Jacob & Pf 2017| |Ruszkowski 1.2017). Reg bly, this body of k i 


29 Just as conduction models probe the fine-grained thermodynamic temperature distribution, the 
coarse-grained mean temperature profile T(x) = fe(x)Te + (1 — fe(x)Th, where fc(x) is the spatially 


varying mass fraction of cold gas, can be reproduced by a mixing length model (Tan et al.|2021 
Chen et al.|2022). 
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simply too enormous to survey here. Instead, we briefly review CR hydrodynamics. Then, 
in keeping with the theme of this section (q3). we consider the specific impact of multi-phase 
gas structure on CR transport. 


3.5.1. CR Hydrodynamics. Since they travel at the speed of light, CRs should zip across 
our Galaxy in ~ 30,000 years. Instead, spallation and radioactive decay products indicate 
CRs have residence times a thousand times longer. CRs originate in discrete, transient 
sources (supernovae) and should be highly directional on the sky. Instead, CRs are remark- 
ably isotropic, to 1 part in ~ 10* (at ~GeV energies, where the energy density peaks). 
These observations make sense if the galaxy is ‘optically thick’ to CRs, so that they scatter 
frequently (^ 10° times before leaving, with ~pc mean free paths) and slowly random walk 
out. Scattering slows the bulk CR propagation speed, generating macroscopic CR pressure 
gradients that push on the gas. 

Magnetic irregularities on scales of order the CR gyroradius (~ 1 AU for a GeV CR) 
can resonantly scatter CRs. These perturbations can originate from two sources: extrinsic 
turbulence, whereby the fluctuations are part of an externally driven cascade, or from the 
CRs themselves (‘self-confinement’), which can amplify Alfvén waves through a resonant 
streaming instability (Kulsrud & Pearce[1969] [Wentzel|1974). While turbulent scattering 
could be important for high energy CRs, theory and observations alike point to a shift 
towards self-confinement for CRs with E « 300GeV (Amato & Blasi| 2018]. The CR 
energy density Ee, which determines CR influence on surroundings, is dominated by low 
energy (~GeV) CRs at the peak of the energy spectrum. For this reason, we focus on the 
physics of self-confined CRs. We emphasize that CR transport is still quite uncertain. For 
instance, there are troubling discrepancies between ‘standard’ models of CR transport and 
observations in the Milky Way [Hopkins et al./2021a), though 
these mostly arises at higher energies. The issue is similar to that with conduction: CR 
scattering rated? ]remain uncertain. Particularly important is the influence of chaotic small- 
scale tangled B-fields, known as ‘Field Line Wandering’ (FLW); see [Mertsch] (2020) for a 
recent review. This renders CR transport effectively diffusive, or even superdiffusive 
Sampson et al.[2022), even if CRs ‘stream’ along field lines. Diffusive CR 
transport can be viewed as a subgrid model for this tangled field structure. We describe 
canonical ‘cosmic ray hydrodynamics’ [Zweibel]2017), but it is on shakier 
ground than say, MHD. 

When tightly self-confined or “coupled”, CRs are trapped by the scattering Alfvén 
waves, which reach an equilibrium amplitude (e.g., B/B ~ 107? for Galactic parameters; 
see between growth and damping processes[!] 
'Their small mean free path renders their collective behavior fluid-like, and their bulk flow is 
described as a sum of advection with the gas, ‘streaming’ relative to the gas with waves at 
the local Alfvén velocity va, and a second order correction from CR diffusioif??] relative to 


the wave frame, due to the finite scattering rate (Skilling]1971). Streaming transport only 


30While these are usually calculated in quasi-linear theory, there has been recent encouraging 


progress in calculating en né rates in Particle-in-Cell and hybrid codes (Bai et al.]2019 2019] [Holcomb] 
& Spitkovsky|2019| [Bai[2022) ). 


EB a SE SS non-linear Landau damping (Cesarsky & Kulsrud| Cesarsky & Kulsrud| Kulsrud}1981), ion-neutral 

damping (De Pontieu et al.|2001), turbulent damping (Farmer & Goldreich|2004), and dust damping 

(Squire et al et al./2021). 
A E this is also referred to as ‘diffusion’, it is distinct from diffusion due to FLW. 
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proceeds down the CR gradient (i.e., to regions with lower CR energy density E.) along 
magnetic field lines. CRs propagating up their gradient damp the streaming instability, and 
do not couple to the gas. 

CRs also push the gas with a force VP. and heat the gas at a rate va - VP. These 
momentum and energy transfer rates can be considerable, but they also depend sensitively 
on local plasma conditions. This last point is important. One reason why CR feedback is 
considered so attractive is that it can act on large scales and suffers weaker losses than say, 
thermal gas. This assumption, which generally comes from low-resolution galaxy-scale simu- 
lations, or when CR streaming is neglected, deserves continued scrutiny in higher-resolution 
calculations. For instance, in the tight-coupling limit, CR pressure traces variations in ve- 


locity, density, and B-field strength (Breitschwerdt et al.|1991): 
P. x (v 4 vA) 4 = (v + B/A/Anp) 4/9, 11. 


Small scale structure in these quantities can thus significantly modulate CR profiles. 
Complicating matters further, CRs are not always coupled to the gas! Coupling is weak 
if wave damping is strong, or if there are insufficient CRs to power the streaming instability. 
CRs also only scatter and couple to the gas if they are anisotropic. As with a radiation 
field, spatial anisotropy implies that CRs with energy density E. have a net flux Fe, with 
some effective drift speed up ~ F./E-. For the streaming instability to be excited, the drift 
speed must exceed the local Alfvén velocity v4. In regions where the CRs are isotropic 
(V.P; = 0), or have small drift speed, vp < va, CRs will not scatter; they decouple from the 
gas and free stream out of these ‘optically thin’ regions at the speed of light. Effectively, CRs 
behave like a radiation field with a very unusual opacity, which depends on the anisotropy 
of the radiation field. This important property, which strongly influences CR transport 
and dynamics, has traditionally been difficult to simulate, as it leads to a rapidly growing 
grid-scale numerical instability, which can only be offset by regularization 
[2010], which is numerically costly (quadratic timestep scaling, At ox (Ax), where Az is 
the grid cell size) . In recent years, two-moment methods adapted from radiative transfer, 
which solve time-dependent equations for both CR pressure P. and flux Fe, have allowed 
for fast, stable and accurate solution of CR hydrodynamics with At œ Ax scaling 


Thomas & Pfrommer Chan et al.|2019). 


3.5.2. CR transport in a multi-phase medium. The CR ‘bottleneck’ is a dramatic example of 
i adi iram or SE Ger Be NA d complexities into D transport 
warm Sking ~ 2 Zi MN gas DRM in m ara ~ » 10°K) gas results in a minimum in 
drift on This produces a ‘bottleneck’ for the CRs: CR density is enhanced as CRs are 
forced to slow down, akin to a traffic jam. Since CRs cannot stream up a gradient, the 
system readjusts to a state where the CR profile is flat up to the minimum in v4; thereafter 
the CR pressure falls again. If there are multiple bottlenecks, this produces a staircase 
structure in the CR profile (Fig [13). A convex hull construction (connecting the highest 
peaks in (v + vA)! with horizontal ridgelines) agrees with numerical simulations 
[et a. [20225). Importantly, since V P. — 0 in the plateaus, CRs there are no longer coupled 
to the gas, and can no longer exert pressure forces or heat the gas. Instead, momentum 
and energy deposition is focused at the CR 'steps'. Small-scale density contrasts can thus 
have global influence on CR driving and heating. Staircase structures appear when CR 
streaming dominates over advection or diffusion. They are seen in 1D simulations of CR 
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Figure 13: The CR bottleneck effect. Left: in steady state, and if CR streaming dominates, 
P, x (va 4- v) œ piel? where the last step assumes v < va and B ~ const. However, 
in the presence of density bumps, this requires CRs stream up their gradient. Right: In- 
stead, CRs ‘bottleneck’ there, and produce a CR staircase structure consisting of plateaus, 
where CRs are decoupled, and sharp gradients, which undergo intense heating and forcing. 
Thus, density fluctuations (e.g., from multi-phase structure) drastically changes the spatial 
footprint of CRs. Figure credit: Navin Tsung. 


driven winds (Quataert et al. 2022a), CR driven acoustic instabilit (Tsung et al.| 2022b), 


and in 2D/3D simulations of thermal instability (Tsung et al.[2022a) and cloud acceleration 
(see below). The importance of bottlenecks depends on the covering fraction of overdense 
gas, and whether they are threaded by B-field lines. 

Cosmic rays can accelerate cold clouds by exerting a direct force. Recall that P.(va + 
v)”: is conserved in the tight-coupling limit (equation [11]. Suppose va > v, so that CR 
streaming dominates. As the Alfven velocity drops within the cloud, CRs bottleneck and 
build up. The resulting steep pressure gradient at the cloud interface accelerates the cloud 
(Wiener otal 20TTa 
et al.12022). Simulations (mostly in 2D) have found efficient acceleration (although it must 
be extrapolated to observed velocities), and that radiative cooling wards off the destructive 
effects of CR. heating. In detail, there are caveats. In warm (T ~ 10*K) clouds, the CRs 
exert their force on the rear end, which is the final bottleneck where they accumulate. 
Thus, they stretch the cloud. Depending on field strength, this differential acceleration 
can shred the cloud (Huang et al.][2022). This mechanism requires field lines to thread 
the cloud, which does happen for clouds which form via thermal instability 
[2022], but does not happen for clouds which undergo magnetic draping. All simulations 
adopt a setup where the hot and cold gas are initially at rest. The bottleneck diminishes 
as their relative velocity v increases, and goes to zero if v is super-Alfvenic (relative to 


33This is an instability in low 8 plasmas where CRs drive sound waves unstable, creating a series 


of closely spaced weak shocks (Begelman & Zweibel|1994). The density jump at the shock serves 


as a propagating bottleneck. 
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the hot gas Alfven speed, va ~ cs,a/v B). Thus, this mechanism alone cannot accelerate 
clouds significantly in a high 8 plasma. Still, even if CRs accelerate the hot gas alone, the 
latter could transfer momentum to the cold gas, producing indirect acceleration. CRs will 
pressurize and broaden the interface; find that the reduction in cooling 
suppresses cool gas growth. 

Self-shielding and modulation of ionized fraction can also impact CR coupling. Ion- 
neutral damping attenuates the MHD waves which couple CR to the gas, and increase 


CR diffusion (Farber et al.| 2018| |Bustard & Zweibel|2021). Also, since va = B/(4mpi)'/”, 


where p; is the ion density, the Alfven speed can rise, rather than fall, in the interiors of self- 
shielded clouds. find that this varying Alfven speed profile creates 
bottlenecks at both the cloud front and back, with modest effects on cloud acceleration 
compared to the single bottleneck case. 


3.5.3. Outlook. This is fast-moving field; our discussion barely scrapes the surface. The 
primary challenge is to settle on the correct equations for CR transport. At present, the 
standard procedure is to consider limits where either CR streaming or diffusion dominates, 
using constant diffusion coefficients calibrated to Galactic constraints. In reality, the lat- 
ter almost certainly vary with environment. Researchers are starting to use models with 
diffusion coefficients calculated from quasi-linear theory; differing assumptions (e.g., for 
wave damping), compounded by non-linear effects, can result in very different outcomes in 


the CGM (Hopkins et al.|2021b). For instance, the low CGM B-fields (and hence Alfven 


speeds) in FIRE simulations mean that CR streaming is unimportant; strong CR diffusion 
must take place to achieve sufficiently rapid CR transport to satisfy y-ray constraints (Chan 
2019), but this may not be true in simulations with stronger B-fields (e.g., 
Voort et al.[2021]. Even if diffusion is simply an effective subgrid model for field line wan- 
dering, how it changes with turbulence and plasma parameters, and the rate of cross-field 
transport, is still unclear. FLW is usually studied in the test particle approximation in a 
static tangled B-field, which ignores the backreaction of CRs (from CR forces and heating) 
on B-field structure, evolution of the B-field due to turbulence, and the constraint that 
self-confined CRs can only propagate down their gradient. Even with the current set of 
equations, unresolved structure (such as multi-phase gas) in simulations can significantly 
change outcomes, as we have seen in our discussion of bottlenecks, and this may also have 
to be handled in a subgrid manner. Finally, increased contact with observations is crucial. 
CRs are currently constrained by non-thermal radio (synchrotron emission from cosmic 
rays electrons) and gamma-ray emission (from pions produced by CR protons); the latter 
more directly constrain CR energy densities. In starburst galaxies, there are well-known 
correlations between radio and gamma-ray luminosities and far-infrared luminosity, which 
probe the star formation rate 
[2012). However, these generally probe conditions at the base of the wind rather than the 
CGM, and constraints on CR transport suffer from degeneracies (e.g., between CR diffusion 
coefficient and halo size; 2011). Comparison with observations requires codes 
with spectral energy resolution (Armillotta et al.|[2021] 
2022). Progress in CR transport could also come from more local, high 
resolution observations (Thomas et al./2020). Given the recent pace of activity in this field, 
there are ample grounds for optimism. 
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4. RESOLVING SMALL SCALES IN LARGE-SCALE SIMULATIONS 


This review has emphasized the multi-scale nature of the multi-phase CGM. The large 
dynamic range on which important physical processes operate means that achieving numer- 
ical convergence in cold gas properties — which generally have smaller characteristic length 
and mass scales than hot gas — can be extremely challenging in large-scale simulations. 
In this section, we compile some key length scales, discuss what is currently achievable in 
cosmological simulations, and indicate potential paths toward progress. When discussing 
convergence, it is important to be specific about the metric of interest. For instance, conver- 
gence in cold gas morphology is much more demanding than convergence in cold gas mass, 
as we explain below. We expect that numerical convergence in column densities, which is 
relevant for observational predictions, is intermediate between convergence in cold gas mass 
and convergence in cold gas morphology, given that it requires convergence in the PDF of 
the 2D projected density field, rather than the full 3D density field of cold gas. In particu- 
lar, for high covering fractions (where a line of sight passes through multiple clouds), small 
scale structure should be washed out in projection. However, exact convergence criteria 
have yet to be quantified (see further discussion in 84.2). 


4.1. Some Characteristic Scales 


We begin by summarizing some key length scales relevant to multi-phase gas in the CGM, 
many of which have been discussed earlier in this article (see also [Li_et_al.|[2020). Since 
many of these are collisional length scales, which scale as A x n^! (for instance, cooling 
lengths | ~ Geet are proportional to Lee) cx n lj, it is useful to define them in terms of 
column density, which is a density invariant quantity. Besides quoting spatial resolution, it 
is useful for simulators to quote resolved cold phase column densities, as this quickly gives 
an idea of which processes can be resolved. We note that while the collisional length scales 
are density invariant, the minimum resolvable columns in simulations in general do depend 
on density (see SES) so this must be taken into account in estimating the range of columns 
resolved in different regions of a simulation. 

Firstly, there are variants of the cooling length l ~ vtcoo1, (where v is a characteristic 
velocity and tcoo1 is a characteristic cooling time), which appear to be necessary to have 
converged cold gas mass (equation[12) or converged cold gas morphology (equation[13). We 
caution that this is an area of active research, and researchers still debate which velocity and 
which cooling time is most appropriate (e.g., see q3.2.1], but some version of these length 
scales is likely to remain relevant when the dust settles. In our opinion, of all the length 
scales listed below, for converged cold gas masses as well as rates of mass, momentum and 
energy transfer between phases, the cold gas survival length scale is the most important 
one to target. 


e Cold gas survival. The survival length is the minimal size for a cloud of overdensity 
x to survive when there is velocity shear relative to the hot gas with Mach number 
M = Glen, Adopting the criterion tcool,mix < tec as in equation (8), so that | ~ 
Usheartcool,mix//X, we obtain: 


—2 Mx100T 3/7 


DE a 108 cm? MT? 12. 


it 18 
Ng ~ 10° cm Hn 
; 


Amix,21.4 
where in the second step we assume isobaric conditions x = pe/pn = Th/Te, A(Tmix) « 
quy (appropriate for 10° < Tmix < 10^? K), and Ta = 10* K. Clouds with Ng > 


mix 
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Ng" can survive and grow in a wind. This criterion is still somewhat uncertain for 
M, x100 >> 1 conditions (X100 = x/100). 

e Sonic length. The sonic length (sometimes referred to as the ‘shattering’ length, see 
93.3.2) is the minimum value of the cooling length cstcoo1, at which the sound crossing 
time tsc ~ l/cs and the cooling time Lee) are equal, so l ~ cstcooi For solar metallicity, 
collisionally ionized cooling curves, it arises at T ~ 10* K, giving 


Ng’ ~ 10" om", 13. 


although it can be larger if cooling is inefficient, e.g. at low metallicity or due to 
photoionization. 


To put these length scales in context, it is useful to compare them against some others. 


For instance, there are length scales associated with the diffusive transport of momentum 
and energy, i.e. viscosity (mediated by ions) and thermal conduction (mediated by elec- 
trons). Viscosity tends to be less important, given long viscous times: the thermal velocity 
of ions is (mp /m.) 


1/2 AO times smaller than the thermal velocity of electrons. Thus, 


simulations which implement Braginskii viscosity in cloud-wind interactions see little effect 
(Li et al.12020). There are two length scales associated with thermal conduction: 


e Electron mean free path. The Coulomb mean free path 
NE? dU cm fmtpTe 14. 


determines the nature of thermal conduction. For Ng > NDA. classical diffusive 
(Spitzer) conduction is applicable. Here, fa: < 1 is a factor which set the reduc- 
tion of the electron mean free path by scattering mechanisms other than Coulomb 
interactions. 

e Field length. The Field length is the length scale at which classical diffusive conduc- 
tion and radiative cooling balance. If the heat flux is given by F = —&(T)VT, then 
V- F n K(D)T/ AL. ~ n? A(T) gives Ap ~ (&(T)T/n? A(T))!?. Conduction (cooling) 
dominates on smaller (larger) length scales. Evaluated in the hot medium, 


Ng ~ 10 7 om^" found Te 5 15. 


where Land < 1 is a suppression factor (e.g., due to tangled B-fields, or additional 


scattering). Static cold clouds with Ny < NE embedded in a hot medium evap- 


orate, while larger ones will condense if pressure exceeds a critical value (McKee &| 


1990). Note that the Field length decreases rapidly with temperature. At 


small scale interfaces between different phases, it sets the local temperature scale 
height, even in the presence of turbulence (Tan & On[2021). Thus, to obtain the cor- 
rect temperature PDF at phase interfaces (e.g., for computing collisional ionization), 
the Field length must be resolved. If not resolved, this can be implemented at the 
subgrid level. 


As of this writing, given their computational expense, we do not believe the electron mean 
free path and the Field length are important length scales to be targeted for typical purposes, 
and if required, they can be incorporated via subgrid prescriptions. 


Finally, there are length scales associated with self-gravity and radiative transfer: 


www.annualreviews.org e Physical Processes in the CGM 


53 


54 


e Jeans length. The Jeans length A; ~ cst is the critical size for clouds to become 
self-gravitating. Evaluated at T ~ 10^ K, it equivalent to a column density: 


Nl ~ 5x 10? cm? PI, 16. 


where P; = n'T/(10? K cm?) = n-2T4 (n-2 = n/(10 ? cm ?)) is the pressure of the 
cold gas. Clouds with Ng > Ni become self-gravitating. Most clouds in the CGM 
are pressure confined rather than self-gravitating. 

e Self-shielding. Most photoionization modeling in the CGM context assumes that 
cold gas is optically thin. However, clouds start to self-shield from hydrogen ioniz- 
ing photons once their column densities start to exceed Nui ~ Gay ~ 1017 cm”, 
where opt is the photoionization cross section at the Lyman edge (13.6 eV). However, 
the total associated hydrogen column density is larger, Nu ~ Niort gr orit where 
THI,crit < 1 is the hydrogen neutral fraction when self-shielding starts to kick in 


(Zheng & Miralda-Escudé|2002). Thus, the column density associated with cooling 
to T ~ 10 — 100K < 10* K, via fine-structure and molecular lines (Robertson A) 
2008} |Krumholz & Gnedin|2011), is larger: 


NES — 10? cm Sa, 17. 


where Zo.1 = Z/(0.1Z). Such large columns are typically only found in the ISM, 
though at high redshift they can be reached also in the inner CGM (e.g. |Stern et al. 
2021b). Thus, most clouds in the CGM will stay at warm (atomic) temperatures 
T ~ 10* K. 


4.2. Achievable Resolutions in Cosmological Simulations 


It is useful to compare the above physical scales for cold gas to the resolution achieved by 
different types of cosmological simulations. We caution upfront that the simple scalings 
below are only crude guides to the scales resolved by different simulations, because this 
depends on the details of the specific simulation and numerical method, including the 
order-of-accuracy of the solver in addition to the cell size or particle mass. As mentioned 
above, which cold gas length scales should be resolved depends on the question asked. For 
example, for many purposes it is likely not necessary to resolve the small sonic/shattering 
length (see discussion in q3.3] and 3.4). 

Most cosmological galaxy formation simulations use a Lagrangian (or quasi-Lagrangian) 
method, in which the mass of resolution is fixed and the spatial resolution adapts according 
to where the mass ends up. This approach has the advantage that it concentrates most of 
the resolution in galaxies, while making it possible to simulate large boxes in which most 
of the volume has low density. A simple estimate of the spatial resolution in the gas for a 
Lagrangian simulation is 


1/3 
ees (=) & 160 pe mè an CAT, 18. 
S , 


where mgas is the mass of gas cells and mgas,3 = Mgas/(10° Mo). In terms of column 
density, the resolution is 


2,1 2 
NES © hgasn £2 5 x 107? cm~? mis una. 19. 
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This shows that lower columns have more stringent resolution requirements (lower mgas) and 
also that the minimum resolved column depends on density in Lagrangian codes. We express 
resolution in terms of a fiducial density n = 107? cm ? because this is representative of cold 
streams in simulations (e.g., Mandelker et al.[2020a), but in reality cold gas has a range of 
densities. The scalings in equations (18) and (19) apply to moving-mesh codes (e.g.,|Springel 

, mesh-free finite mass (MFM) or mesh-free finite volume (MFM) methods (Hopkins| 
2015), and adaptive mesh refinement (AMR) codes using quasi-Lagrangian refinement (e.g., 


od Geen ean SPH codes are also widely used 


for cosmological simulations but the spatial resolution for these is better approximated 


by the size of the smoothing kernel, which typically contains Nnpg ~ 32 — 100 particles. 
Adopting a fiducial neighbor number Nnbg œ~ 60 (e.g., oe the spatial 
resolution of an SPH simulation can be estimated as ~ Ni 4x coarser than implied by 
equation (18), with the cell mass replaced by the SPH Be mass. The minimum resolved 
column densities are correspondingly larger. The fiducial cell mass mgas © 1,000 Mo in 
the equations above are representative of today’s state-of-the-art zoom-in simulations for 
Milky Way-mass galaxies (e.g., [Grand et al.|2021). However, full cosmological boxes with 
side lengths ~ 100 Mpc typically have a much larger cell mass Mgas ~ 10° Mo (e.g., 


[et a [2014] Springel et al. Springel et al 2018), corresponding to a ~ 10x coarser spatial resolution, i.e. 


~ 1.6 kpc or Ng ~ 5 x 10/9 cm ? in n = 107? cm7 3 gas. 


Comparing the above resolution scalings with characteristic scales in q4.1|suggests that 
state-of-the-art zoom-in simulations can resolve the Jeans scale of cold gas (eq. and 
may possibly start to marginally resolve the cold gas survival scale (eq. |12), but the 
sonic/shattering scale is generally not resolved. This implies that the highest resolution 
zoom-in simulations are approaching resolutions where the cold gas masses are converged, 
but that the cold gas morphology (e.g., the amount of small-scale structure) remains un- 
converged. This is consistent with the moving-mesh, zoom-in simulations of 
[et al.| (2019), who find that the total HI mass in the CGM is converged across the resolution 
levels investigated, but that the amount of small-scale structure in HI maps increases with 
increasing CGM resolution. Consistent with the varying small-scale structure, these authors 
find that the covering factors of HI above different column density thresholds also depend 
on resolution. Large-volume cosmological simulations resolve the Jeans scale but neither 
the cold gas survival scale nor the sonic scale. Thus, cold gas masses and morphology may 
both be unconverged in such simulations, although the degree of convergence likely depends 
on subgrid models and whether they are “recalibrated” at different resolutions (e.g., 
ot al [2015). 

In the last few years, multiple groups have noted that a Lagrangian approach may 
not be ideal for CGM studies because the concentration of resolution in galaxies implies 
a relatively low spatial resolution in dilute gas outside galaxies. To address this, new 
refinement schemes have been devised to enforce either a minimum or an uniform spatial 
resolution in a specified region EE into the CGM or E 'This approach has been 
pe in ad 2019] 2019} [Peeples et al.|2019b) and moving-mesh codes 
2019), and it Pate in principle also be used in mesh-free codes that 
Mmm splitting of T. elements. In these simulations, the spatial resolution Az is 
fixed while the mass resolution depends on density as 


Mres = (Axr)*p ~ 2.5 x 10° Mo (Ax) }n_2 20. 


www.annualreviews.org e Physical Processes in the CGM 


55 


56 


and the minimum resolved column density is 
Ng? = (Aa)n 2: 3 x 10? cm? (Ax)in_2, 21. 


where (Az); = Az/(1l kpc). In published simulations with uniform or enhanced CGM 
refinement, the maximum cell size is (Az); ~% 0.5 — 1 (quoted here at z = 0, since some 
simulations prescribed the maximum cell size in proper units and others in comoving units). 
Therefore, for the fiducial cold gas density n = 107? cm^?, the minimum resolved column 
is comparable to (or larger than) what is achievable with Lagrangian zoom-ins (eq. [19]. 
These studies focused on CGM resolution have been influential in demonstrating clearly 
that finer structures are resolved as the spatial resolution is increased, especially in the 
cold gas (see an example of a resolution study in Fig. E. However, in agreement with 
the analytic estimates here, the cold gas morphology is not converged even at the highest 
spatial resolutions simulated to date. An important limitation of simulations that limit the 
size of resolution elements to a prescribed maximum is that while this produces very fine 
mass resolution in low-density gas, this tends to be predominantly in hot gas rather than 
in the cold, dense gas. Unfortunately, it is not clear how much higher resolution in hot 
gas helps, relative to Lagrangian refinement schemes that achieve comparable resolution in 
cold gas but coarser resolution in hot gas, since the hot gas tends to be much smoother. A 
promising variant may be to adopt refinement criteria that are tailored to specific processes 
that one wishes to resolve, such as refining on the local cooling length (e.g., 
[et al.[2021). The insights on relevant physical scales from small-scale studies will play an 
important role in designing such refinement criteria. 


4.3. Incorporating Small Scales via Subgrid Recipes 


Another approach is to distill the insights from small-scale simulations as subgrid models. 
There has been some progress on this front for cloud-wind interactions. For instance, PhEW 
(Physically Evolved Winds; calculates the exchange of mass and met- 
als between an SPH cloud ‘particle’ and its surroundings, incorporating isotropic thermal 
conduction, hydrodynamic instabilities, and the compressive and elongational effects of a 
bow shock, calibrated to the high resolution cloud crushing simulations of[Brüggen & Scan-| 
(2016), albeit only in the cloud destruction regime. PhEW has been implemented 


in cosmological simulations (Huang et al.|2022). It improves numerical convergence, since 


small-scale interactions are handled in subgrid. 
While omitting some of the physics in PhEW, such as conduction and the effects of a bow 


source terms for cloud entrainment and growth (as discussed in 83.2.1), which reproduce 
the results of idealized wind tunnel and shearing layer simulations in this regimd?*| 
include a wind thermalization term (see also 2021) 
which matches the rising entropy seen in high-resolution wind simulations (Schneider et al.| 
2020). Their model enables fast exploration of parameter space in semi-analytic models of 
galaxy formation, and there are plans to incorporate it into hydrodynamic simulations. 

A related approach is to use a two-fluid approximation, similar to that used in cosmic 
ray or radiation hydrodynamics simulations, where hot and cold gas exchange mass, energy 


34 As noted in d3.2.4] cloud survival and growth in the high overdensity or Mach number regimes 
is still not completely clear, so their extrapolations in these regimes is necessarily uncertain. 
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and momentum via prescribed source terms (Weinberger & Hernquist|2022). Of course, 


all the magic lies in these source terms, which must encode our growing understanding of 
multi-phase interactions. These prescriptions have yet to be formulated in this approach. 
These developments illustrate how small and large scale simulations can have a fruitful 
synergy, similar to treatments of the ISM and star formation in zoom-in or cosmological scale 
simulations. While this approach holds considerable promise, these are early days. More 
work is needed, both in understanding the physics of small-scale interactions and encoding 
these insights into subgrid recipes. Perhaps the mark of a mature subgrid recipe would 
be the ability to match high-resolution meso-scale simulations in disparate environments, 
such as multi-phase galactic winds (Schneider et al.|2020) or ‘chaotic cold accretion’ onto 


black holes (Gaspari et al.2013 2018). Such a model could then be used with confidence 


in large-scale simulations. 


5. OUTLOOK AND FUTURE DIRECTIONS 


While galactic halo gas has been studied for decades, the advent of the HST Cosmic Origins 
Spectrograph in the late 2000s transformed the field. Many researchers (present authors 
included) voted with their feet and stampeded to work on the CGM. [Tumlinson et al.|(2017) 
reviewed the remarkable ensuing observational progress in this journal. Their Figure 1, 
featured in innumerable conference talks, is a wildly influential cartoon of gas flows in the 
CGM, showing filamentary accretion from the IGM, bipolar winds, and gas recycling. In 
their memorable formulation, "he CGM] is potentially the gas fuel tank, waste dump, 
and recyling center all at the same time." Gas flows indeed lie at the heart of the CGM's 
influence on galaxy formation. The review concluded that the ‘missing’ baryons/metal 
problem, which had plagued researchers for years, could largely be laid to rest. 

'This review partly serves as a theoretical counterpart, summarizing theoretical progress 
on understanding gas flows — both on halo scales, as well as small-scale mass transfer 
between phases. While the [Tumlinson et al cartoon picture is still informative, the 
correct physical picture for the CGM likely evolves with both redshift and halo mass, with 
additional processes contributing. For instance, besides accretion via cold streams, we have 
reviewed the potential impact of hot mode accretion, thermal instability and precipitation, 
cooling-induced condensation onto cold gas ‘seeds’, and the contribution of ‘intergalactic 
gas transfer’ between galaxies. Our new cartoon picture in Figure[I]illustrates some of this 
complexity, and in particular emphasizes the multiscale structure expected in cold gas. 

What are potential agenda items for continued theoretical progress? They include: 


e Clarity on cold vs. hot mode accretion. Despite more than two decades of work, the 
role of cold vs. hot mode accretion is not yet fully settled. A key issue we have 
reviewed is progress on the survival of cold gas streams and clouds, processes which 
are not well resolved in existing cosmological simulations, using idealized simulations. 
Another issue is the effects of different accretion modes on the properties and evolution 
of central galaxies. While there has been much attention devoted to the role of cold 
streams in feeding high-redshift disk galaxies (which are highly turbulent and may be 
considered “thick” ), recent work suggests that hot-mode accretion may be important 
for the formation of thin disks (see (2.1.6). Future work should aim to clarify how the 
mode of gas accretion affects galaxy formation, and how this depends on mass and 
redshift. One specific area where a more systematic analysis of simulations would be 
beneficial would be the mechanisms that exchange angular momentum in the CGM, 


www.annualreviews.org e Physical Processes in the CGM 


58 


and how these exchanges affect the delivery of angular momentum to galaxies (82.2). 
Understanding the structure of cold gas. Most work on small-scale structure in cold 
gas is based on high-resolution simulations with idealized initial and boundary con- 
ditions. What are the properties of the cold gas in realistic galaxy halos? Are halos 
filled with a quasi-uniform cold ‘mist’, or composed of discrete cloud complexes, as 
seen in HVCs? If there is structure on all scales, what is the physics that determines 
the cloud mass function, including its low and high mass cut-offs? And what mass 
range is most germane to particular observational probes? 

Connecting large and small scales. How do we build on recent theoretical progress 
on large and small scales to develop models of the CGM that realistically incorporate 
the cosmological environment, and still have enough physics/resolution to model fine- 
scale structure in the cold gas? Converged cold gas masses and column densities 
at CGM scales is one goal. Another would be models that can be meaningfully 
compared with high-resolution observations of kinematic profiles, which are sensitive 
to detailed structure along the line of sight. Arguably, neither existing cosmological 
nor idealized simulations simultaneously have all the physics and resolution to do this 
well at present. 

Progress om non-thermal physics. B-field strengths affect MHD forces, turbulence, 
and CR streaming speeds. The MHD equations are in principle straightforward to 
solve. However, in practice, B-fields on CGM scales vary greatly between different 
simulation groups, likely because magnetic amplification and advection depend on 
the feedback model. Subgrid models for unresolved tangling might be needed, given 
its effects on CR transport and thermal conduction. Unlike MHD, it is as yet unclear 
whether we are even solving the correct equations for CR transport. Progress will 
likely come with continued confrontation of observations with physical models. 


What observational inputs would most benefit theorists? A partial inventory includes: 


e Better characterization of the hot CGM. In most galaxies at L, and below, X-ray emis- 


sion from coronal gas is too faint to be observed, and we are confined to observing cold 
gas (T ~ 104 K). Even if pressure balance holds between phases, the unknown contri- 
bution of non-thermal pressure support and large uncertainties in cold gas properties 
inferred from ionization model modeling (e.g., 
imply considerable uncertainties in hot gas pressures. Direct constraints on coronal 
gas properties — which affect mass, momentum and energy exchange between phases 
— and their radial profiles are sorely needed. FRB dispersion measure measurements 
(which probe the total columns of free electrons), as well as a combination of thermal 
and kinetic Sunyaev-Zeldovich measurements (which together can reveal density and 
temperature profiles; e.g.,|Schaan et al. 2021), could potentially shed light on this. 

Better characterization of non-thermal components. We have reviewed the significant 
impact that B-fields, cosmic rays, and turbulence play on physical processes in the 
CGM. At the same time, depending on underlying assumptions, simulations can give 
widely disparate predictions for non-thermal energy content, from energetically in- 


significant to magnetically or CR dominated halos (e.g., Hopkins et al.|2020 


2020). FRBs could shed more light on B-fields and turbulence, as well as small-scale 


clumping (Prochaska et al.|2019| ?). The CR component may be the most difficult 


— while in principle it is constrained by gamma-ray or radio observations, in practice 
these are only detectable near the host galaxy, rather than far out in the CGM. 
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e Spatially resolved observations. Most information about the CGM comes from quasar 
absorption line data which suffers from poor spatial coverage (typically just one pencil 
beam per halo); what is needed is data that combines high spatial sampling and 


high spectral resolution (Peeples et al.|20192). Gravitational lensing has been used 
fruitfully to enhance the density of background quasars (e. e Chen et al.[2014] [Bowen] 
[Lopez et al.[2018], and future studies using 30m-class telescopes will take 
spectra of multiple fainter background galaxies per halo (Steidel et al.|2010). Even 


better, emission lines could provide detailed spatial and kinematic maps. At high 
redshift, halos are already seen glowing in Lyo emission (82.1.5) but the complexity 
of Lya radiative transfer precludes simple interpretation. Optically thin metal lines, 


although faint, are potentially much more informative 
[Piacitelli et al.|2022]. Directly imaging the CGM was highlighted as a 
key discovery area in the Astro 2020 decadal survey 
[& Medicine[2021). Although CGM imaging has recently been enabled by new integral 


field spectrographs on 8-10m telescopes, imaging diffuse gas is presently limited either 
to the brightest, inner regions of the CGM, or to special cases where the gas is 
illuminated by a quasar. Order-of-magnitude improved sensitivities will be needed to 
map the CGM of normal galaxies on larger scales. 


A potential preview of the future of CGM studies comes from their more massive sib- 
ling, the ICM. All the same physical ingredients are present: radiative cooling, turbulence, 
stratification, multi-phase gas, B-fields and cosmic-rays. However, we see the hot medium 
directly, in spatially resolved X-ray observations. We also see cold filaments directly in 
resolved Ha and CO emission. Despite short central cooling times, cooling flow models 
(Fabian|1994) were shown by Chandra and XMM observations to be untenable in clusters 
(Peterson & Fabian|/2006). Instead, the ICM is understood to be in approximate hydro- 
static and thermal equilibrium — not too dissimilar from a star, with feedback from the 
central AGN playing the role of nuclear reactions (McNamara & Nulsen|2007). We see the 
feedback directly, in the form of spatially resolved AGN blown bubbles. Non-thermal com- 
ponents are much better characterized: ~ uG magnetic fields (leading to 8 ~ 100 plasmas) 
are constrained by Faraday rotation; cosmic rays (P./P X few%) by lack of gamma-ray 
emissio] and kinematics of highly subsonic M ~ 0.1 — 0.2 turbulence by emission line 
broadening (Hitomi Collaboration et al.|2016). Bulk sloshing motions are constrained by 
observed eene e Eee pun] Which is not to say that matters 
are fully resolved. For instance, theorists still debate how AGN energy is transmitted and 
isotropized in the ICM; conduction and viscosity, and the role of kinetic instabilities in 
modulating them, are highly uncertain; the origin of cold filaments, as well as their kine- 
matic and energetic coupling with the hot phase, is still debated (Donahue & Voit|2022). 
But the field is certainly in a significantly more mature state than CGM studies, due to the 
transformative power of spatially resolved observations. 

One possible lesson for CGM theorists is, wherever possible, to use systems with a 
wealth of high resolution data as laboratories for models of key physical processes, before 
sallying forth into terra incognita. If, after multiple decades there is no direct evidence for 
thermal conduction in the ISM or ICM, it is hard to imagine we could detect its effects in 


35This was considerably smaller than expected (Pinzke & Pfrommer|2010), given the expected 


high efficiency of CR acceleration at strong accretion shocks, and the fact that the ICM is largely 
a closed box — a sobering reminder of surprises and uncertainties in CR physics. 
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the CGM anytime soon. A powerful laboratory for physical models of CR transport is our 
Galaxy Hopkins et al.12021a). Models of small scale structure 
in photoionized gas should confront ~ 1 — 10^ AU scale atomic structures and extreme 
scattering events in the ISM (Stanimirovié & Zweibel|2018), which may also constrain MHD 
turbulence on CR gyro-scales. Models of turbulent mixing layers and cold gas condensation 
onto cometary tails must confront resolved observations of HVCs and cluster filaments, 
jellyfish galaxies (Müller et al.[2021), and AGB winds such as Mira A 
Li et al.[2019). 

It’s also good to remember that progress is often fueled by cross-fertilization from other 
fields. Precision cosmology was greatly aided by the adoption of Markov Chain Monte 
Carlo, and machine learning techniques are playing a similar role today 
[2022). Similarly, the scale problem that CGM theorists confront is not too dissimilar from 
the problem of thermo-nuclear burning on white dwarfs (Schmidt et a1. [2006], which is even 
more difficult to treat due to the high temperature sensitivity of nuclear burning rates. It 
is also similar to what engineers working on turbulent mixing face (Sreenivasan] 2019). 
Climate science simulations tackle another complex multi-scale problem: indeed, clouds 
are the largest source of uncertainty and usually have to be handled in a subgrid fashion 
(Huang et al.|2022). Even if the stakes for humanity are not quite as high, our fledging 
efforts to model galactic atmospheres can draw hope and inspiration from the progress in 
climate science simulations over the past few decades. 
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ACRONYMNS AND DEFINITIONS 


Virialized gas: Gas with thermal velocities set by the gravitational potential. In massive 
halos, this corresponds to a hot phase (Tvir Z 10° K). 


Cold flows/cold streams: Cold halo gas, typically at a temperature T ~ 10* K and 
filamentary in morphology. Associated with “cold mode” accretion. 


Cooling flow: Model for hot gas accretion in halos driven by radiative cooling and 
neglecting feedback. Associated with “hot mode” accretion. 


Precipitation: Formation of cold gas via linear thermal instability in a gravitationally 
stratified hot medium. The overdense cold gas eventually falls, like rain. 


Condensation: Conversion of hot to cold gas when initially hot gas mixes with cold gas, 
and the mixed gas cools. The hot gas ‘condenses’ onto cold gas seeds. 


Turbulent mixing layer: A turbulent interface between hot and cold phases where 


the interaction between turbulence and radiative cooling governs the exchange of mass, 
momentum and energy between phases. 
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